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ABSTRACT 

Context. The disks around Herbig Ae/Be stars are commonly divided into group I and group II based on their far-infrared spectral 
energy distribution, and the common interpretation for that is flared and flat disks. Our understanding of the evolution of these disks 
is rapidly changing. Recent observations suggest that many flaring disks have gaps, whereas flat disks are thought to be gapless. 
Aims. The different groups of objects can be expected to have different structural signatures in high-angular-resolution data, related 
to gaps, dust settling, and flaring. We aim to use such data to gain new insight into disk structure and evolution. 

Methods. Over the past 10 years, the MIDI instrument on the Very Large Telescope Interferometer has collected observations of 
several tens of protoplanetary disks. We modeled the large set of observations with simple geometric models and compared the 
characteristic sizes among the different objects. A population of radiative-transfer models was synthesized for interpreting the mid- 
infrared signatures. 

Results. Objects with similar luminosities show very different disk sizes in the mid-infrared. This may point to an intrinsic diversity 
or could also hint at different evolutionary stages of the disks. Restricting this to the young objects of intermediate mass, we confirm 
that most group I disks are in agreement with being transitional (i.e., they have gaps). We find that several group II objects have 
mid-infrared sizes and colors that overlap with sources classified as group I, transition disks. This suggests that these sources have 
gaps, which has been demonstrated for a subset of them. This may point to an intermediate population between gapless and transition 
disks. 

Conclusions. Flat disks with gaps are most likely descendants of flat disks without gaps. Potentially related to the formation of 
massive bodies, gaps may therefore even develop in disks in a far stage of grain growth and settling. The evolutionary implications 
of this new population could be twofold. Either gapped flat disks form a separate population of evolved disks or some of them may 
evolve further into flaring disks with large gaps. The latter transformation may be governed by the interaction with a massive planet, 
carving a large gap and dynamically exciting the grain population in the disk. 

Key words, protoplanetary disks, techniques: interferometric, planet-disk interactions, stars: pre-main sequence 


1. Introduction 

Detailed study of the structure and evolution of protoplanetary 
disks is a necessary prerequisite for understanding planet forma¬ 
tion. The dust- and gas-rich protoplanetary disks set the bound¬ 
ary conditions for the formation, initial dynamics, composition, 
and even the actual presence of planetary systems. Planetary sys¬ 
tems are observed to be extremely diverse, so that understanding 
this diversity translates into constraining the onset of planet for¬ 
mation within the disks. 

Objects that have been under special scrutiny are the 
intermediate-mass young objects known as Herbig Ae/Be stars 
(Herbig 1960[ Waters & Waelkens 1998| l. From the point of view 
of star formation, they present the link between the lower-mass 
T Tauri objects and the massive young stellar objects. From the 
point of view of planet formation, on the other hand, they rep¬ 
resent the progenitors of debris-disk objects around A-type stars 
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with or without detected pl anetary systems, such as Vega, j3 Pic, 
Fomalhaut, and HR4796A ( jHolland et al.|l998| n 

Our understanding of the evolution of Herbig Ae/Be star 
disks (or, Herbig Ae/Be disks), itself, is in evolution. The idea 
that the disks are gradually dispersed and the central objects 
evolve into “naked” main-sequence stars is well established 
(e.g., |Lada |1987| l. Additionally, the potential mechanisms that 
drive the dispersion process have probably been identified (see, 
e.g., the overview in Alexander et al.|2014| l. However, the cou¬ 
pling of these dissipation mechanisms to observational signa¬ 
tures for the global structure of the disks is highly non-trivial. 


* It is interesting to note that several of the few directly imaged plan¬ 
etary companions are found around A-type stars, which are descen¬ 
dants of Herbig Ae/Be stars. Examples of A-type exoplanet host stars 
are HR 8799 (e.g., |Marois et al.|2008^ , HD 9508 6 (e.g., [Rameau et ah] 
|2013j , A-And (e.g., [Carson et al.|2013^,/8 Pic (e.g., [Lagrange et al.|2010^ , 
and HD 100546 (e.g., Quanz et al. 2013^. 
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Herbig Ae/Be disks 




Waelkens et al. (1994) 


Meeus et al. (2001) 


Maaskant et al. (2013) 


Fig. 1. Different evolution scenarios proposed for Herbig Ae/Be disks. The Herbig Ae/Be phase is indicated in the red boxes. 


We give an overview of the proposed evolutionary scenarios for 
Herbig Ae/Be disks in Fig. and describe them below. 

Waelkens et al.| (|1994|l have distinguished a class of Herbig 


Ae/Be stars with a broad dip in the 10-//m region of their infrared 
excess. This double-peaked spectral energy distribution (SED) 
was interpreted as representing a transition phase between clas¬ 
sical broad infrared excesses of young stars and faint cool ex¬ 
cesses of debris-disk objects. The conjectured explanation for 
this dip by these authors was the development of a physical gap 
in the radial dust distribution, which is also the interpretation of 
n den Ancker et^TI l 1997| l. The origin of this gap would be the 
formation of larger bodies within the disk. The full evolutionary 
scenario that then emerged was shown in Malfait et al. ( 1998| l. 

A new classification of the Herbig Ae/Be objects has been 
proposed by |Meeus et al.| ( |2b01) l. Based on the shape of its mid- 
to far-infrared continuum, a source was classihed as “group 11” 
when the continuum was like a power law and “group 1” when 
an additional cold component was present on the power-law con¬ 
tinuum. The main hypothesis for the spectral behavior is a mor¬ 
phological difference: group 1 sources having flaring disks, and 
group 11 sources have flat (or self-shadowed) disks. The physi¬ 
cal origin for this difference could be grain growth (Dullemond 
& Dominik|2004a] l, decreasing the opacity throughout the disk, 
and/or grain settling ( [Dullemond & Dominik|2004b) l, decreasing 
the irradiated surface of the disk. A likely evolutionary scenario 
for disks, where grains are expected to be growing and settling, 
is that group 1 sources evolve into group 11 sources. 

These evolutionary scenarios are based entirely on spatially 
unresolved observations. High-angular-resolution mid-infrared 
imaging of Herbig Ae/Be objects is challenging the observa¬ 
tional picture. Honda et al. p012| l point out that many group 
1 sources are found to have a gapped disk. A new evolutionary 
scenario proposed by [Maaskant et al.] ( |2013| l is that group 1 and 
group 11 sources are different successors of a common ances¬ 
tor: a primordial flaring disk. Gap formation would then have 
preceded the collapse of the outer disk in transitional group 1 
sources, whereas in group 11 sources, grain growth and settling 
would have flattened the outer disk. 

The notion that all group 1 sources may have gaps confirms 
the importance of spatially resolved observations. Typical radii 
of directly detected gaps are > 20au (Maaskant et al. 2013| l, 
which is close to the resolution limit of the observations (for a 
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typical distance of 150 pc). Moreover, for the group 11 sources, 
the amount of structural information from direct observations is 
very limited, since the disks tend to appear faint and small. An 
interesting alternative to direct imaging is mid-infrared interfer¬ 
ometry: 

1. the angular resolution for typical observations is a factor of 
10 higher than for direct imaging; 

2. the mid-infrared wavelength range corresponds to the ther¬ 
mal emission of small dust grains in the inner 1 - 10 au (scal¬ 
ing with the stellar luminosity, see [Henning & Meeus[2011) , 
which is a good tracer for the disk structure; 

3. spectrally resolved observations in the 8 - 13 micron atmo¬ 
spheric window provide additional information on the com¬ 
position of the disk’s small grain population. 

Mid-infrared interferometry was successfully used to resolve the 
gaps in the group 1 sources HD 100546 (~ 13 au;[Benisty et al.[ 
[ 20 T 0 I 1 , HD 139614 (~ 6 au; [Matter et al.|2014[ ), and HD 179218 
(~ 10 au; Menu et al., in preparation). 

Given this new evolutionary scenario for Herbig Ae/Be 
disks, many new questions need to be addressed. Are group 1 and 
11 sources really two distinct classes? If a “common ancestor” for 
the two groups exists, what is its nature? Do both groups evolve 
into classical debris-disk objects and/or planetary systems? 

In this work, we aim at investigating the structural differ¬ 
ences within Herbig Ae/Be disks, and protoplanetary disks in 
general, from a large set of mid-infrared interferometric data. 
During its ten years of operation, the MlD-infrared Interfero¬ 
metric instrument (MIDI) on the Very Large Telescope Interfer¬ 
ometer (VLTl) has been used to observe several dozens of proto¬ 
planetary disk objects. This work presents a statistically relevant 
compilation of targets observed with this instrument. In Sect. 
we give an overview of the sample, the observations, and the data 
reduction. Section [^discusses the model choice for interpreting 
these data and the modeling results are presented in Sect. In 
Sect. 1^ we present a size-luminosity relation for the full set of 
protoplanetary disks. This relation covers a broad range in stel¬ 
lar properties, which is interesting but also possibly limits the 
diagnostic power for structural differences within a subgroup. 
We limit the further analysis in Sect. |^to the Herbig Ae stars in 
the sample. Finally, we discuss the results for this group of stars 
in Sect. [ 7 ] in terms of the evolution of the disks. A summary and 
conclusions of this work can be found in Sect.jSj 
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Table 1. Overview of the sample and relevant stellar properties. HAe/HBe = Herbig Ae/Be star, TT = T Tauri star, embHAe = embedded Herbig 
Ae star, DD = debris disk object. 


# 

name 

R.A. (J2000) 

Dec. (J2000) 

type 

d 

7’eff 

L 

references 



(h m s) 

(O , 


(pc) 

(K) 

(Ls) 


1 

T,kHo-330 

034548.28 

+32 2411.9 

TT 

250 ± 50 

5800 

11/5 

1 , 2 ,a 

2 

V892Tau 

0418 40.61 

+28 1915.5 

HAe 

142 ± 14 

11220 


3,b 

3 

RYTau 

04 2157.41 

+28 2635.5 

TT/HAe 

142 ± 14 

6310 

15/« 

3,b 

4 

LkCal5 

04 3917.79 

+22 2103.4 

TT 

142 ± 14 

4350 

0.8+^-2 

3,c 

5 

DRTau 

0447 06.20 

+ 165842.8 

TT 

142 + 14 

4060 

3,c 

6 

GM Aur 

04 55 10.98 

+302159.5 

TT 

142 ± 14 

4730 

1 

-Q/ 

60+1 

11+T 

23^1 

49+1 

3,c 

7 

AB Aur 

04 55 45.84 

+30 33 04.2 

HAe 

139 ± 19 

9800 

4,d 

8 

SUAur 

04 55 59.38 

+30 3401.5 

TT/HAe 

146 ± 50 

5860 

4,c 

9 

HD 31648 

04 58 46.26 

+295036.9 

HAe 

137 ± 25 

8200 

4,e 

10 

UXOri 

05 0429.98 

-03 47 14.2 

HAe 

460 ± 50 

8710 

5,6,f 

11 

HD 36112 

05 3027.52 

+25 19 57.0 

HAe 

279 + 70 

7800 

4,e 

12 

HD 36917 

05 3446.98 

-05 34 14.5 

HAe 

375 + 30 

10000 

230""° 

7,e 

13 

CQ Tau 

05 35 58.46 

+244454.0 

HAe 

113 ±23 

6750 

5 4+^9 

-2,-9 

23+3 

69+1 

4,b 

14 

V1247 0ri 

05 38 05.24 

-01 1521.6 

HAe 

385 ± 15 

7250 

8,f 

15 

HD 38120 

05 43 11.89 

-045949.8 

HAe 

375 + 30 

11000 

7,e 

16 

/?Pic 

0547 17.08 

-5103 59.4 

DD 

19.44 ± 0.05 

8052 

8^« 

4,g 

17 

HD 45677 

0628 17.42 

-1303 11.1 

HBe 

279 + 73 

21380 

600+;*“ 

4,h 

18 

VYMon 

063106.92 

+ 102604.9 

embHAe 

800 + 300 

12023 

3800++^!; 

9,10,1 

19 

HD 259431 

0633 05.19 

+ 101919.9 

HAe 

660 + 100 

14000 

760+1f 

7,e 

20 

21 

RMon 

HD 50138 

063909.94 

065133.39 

+08 4409.7 
-0657 59.4 

HBe 

HBe 

760 ± 300 
392 + 86 

30903 

15490 

2500/|5« 

2500+}™ 

io.j 

4,h 

22 

HD 72106 

08 29 34.89 

-38 3621.1 

HAe 

279 + 88 

8750 

42+^^ 

4,k 

23 

HD 87643 

1004 30.28 

-58 39 52.0 

HBe 

1500 + 500 

17000 

41000+f™ 

11,1 

24 

CRCha 

105906.97 

-77 0140.3 

TT 

160+ 15 

4900 

3.0+;}“ 

12,m 

25 

HD 95 881 

110157.61 

-713048.3 

HAe 

170 ±30 

8990 

^ 11 

13,n 

26 

DICha 

11 07 20.72 

-77 3807.2 

TT 

160+ 15 

5860 

9.5+f} 

12,m 

27 

HD 97048 

1108 03.32 

-77 3917.4 

HAe 

158+ 16 

10000 

40+}^ 

4,h 

28 

HPCha 

1108 15.09 

-77 33 53.1 

TT 

160 ± 15 

4205 

T 9+^ 5 

1 5 

12,m 

29 

FMCha 

110953.40 

-76 3425.5 

TT 

160 ± 15 

4350 

5.8+‘-‘ 

12,m 

30 

WWCha 

11 1000.10 

-763457.8 

TT 

160 ± 15 

4350 

6 5+^-’^ 

u.D 2 3 

12,m 

31 

CV Cha 

11 12 27.70 

-764422.3 

TT 

160 ± 15 

5410 

4 6+*-^ 

H-.U , ^ 

12,m 

32 

HD 98922 

112231.67 

-53 2211.4 

HAe 

1150±515 

10500 

5400+5"*™ 
J+uu „oo 

14+j 

4,e 

33 

HD 100453 

11 33 05.57 

-5419 28.5 

HAe 

122 ± 10 

7400 

4,n 

34 

HD 100546 

1133 25.44 

-701141.2 

HAe 

97 ±4 

10500 

24+5 

1 .2/M 

4,n 

35 

TCha 

115713.55 

-7921 31.5 

TT 

108 ±9 

5890 

14,h 

36 

HD 104237 

120005.08 

-7811 34.5 

HAe 

115±5 

8410 

49+5 

4,n 

37 

HD 109085 

12 32 04.22 

-161145.6 

DD 

18.28 ±0.06 

6784 

5-1-02 

4,g 

38 

DKCha 

12 53 17.23 

-77 0710.7 

embHAe 

178 ± 18 

7200 

19+} 

12,0 

39 

HD 135344 B 

15 15 48.43 

-37 0916.0 

HAe 

142 ± 27 

6750 

12+1 

15,e 

40 

HD 139614 

154046.38 

-42 29 53.5 

HAe 

142 ± 27 

7600 

9 - 3/37 

15,d 

41 

HD 142666 

15 5640.02 

-22 0140.0 

HAe 

140 ± 20 

7900 

19/5 

7,16,e 


2. MIDI sample 


2.1. MIDI and the interest of mid-infrared interferometry 


In 2002, MIDI ( |Leinert et al.|2003] l was installed at the VLTI in 
Chile, and it became a unique instrument for its combination of 
high spatial resolution and spectral resolution {R - 30 - 230) in 
the N band (T = 8 - 13 ;um). For protoplanetary disks, MIDI is 
sensitive to the emission of small dust grains in the inner 1 -10 au 
of the disk (assuming a typical distance of 150 pc). For various 
reasons, this region is highly interesting and contains the imprint 
of ongoing physical processes. With MIDI, it has been shown 
that characteristic sizes of disks can be matched with the inner- 
disk geometry ( [Leinert et al.||2O04) l. Another important hnding 
is that the inner parts of protoplanetary disks can be highly crys¬ 
talline ( |van Boekel et al.| 2004[ ). Finally, the MIDI-detected ge¬ 
ometry of the inner rim of disks can possibly be linked to the in¬ 


teraction with a planetary companion (Mulders et al. 2013 Menu 
let al.|2014| l. 

These discoveries show that MIDI provides a useful and nec¬ 
essary counterpart for high-spatial-resolution observations that 


are also sensitive to large grains (e.g., ALMA) or much hotter 
dust (e.g., PIONIER, AMBER). 


2.2. Sampie and data overview 

Over the past ten years, over 100 young stellar objects with disks 
have been observed with MIDI. In this work, we focus on the 
intermediate-mass objects with, in general, low optical extinc¬ 
tion. The objects are a compilation of sources from longer lists 


of (candidate) Herbig Ae/Be stars in 

The et al.|(|1004|l,|Malfaiti 

let al.| (|19981l, and |Vieira et al.| (|2003 

1. In total, data sets of 38 


Herbig Ae/Be stars are included in the sample. 


The Herbig Ae/Be stars represent only part of the young stel¬ 
lar objects (YSOs) that are observed with MIDI. jBoley etak] 
( |2013| l present an overview and analysis of a large sample of 
massive YSOs observed with MIDI. The transition between the 
high-mass end of the Herbig Ae/Be stars and the genuine mas¬ 
sive YSOs is not clear-cut, and we include a handful of probably 
more massive targets in the sample. 
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Table 1. continued. 


# 

name 

R.A. (J2000) 

Dec. (J2000) 

type 

d 

7’eff 

L 

references 



(h m s) 

(O , 


(pc) 

(K) 

(io) 


42 

HD 142527 

155641.88 

-42 19 23.2 

HAe 

233 + 53 

6260 

50+^^ 

4,n 

43 

HD 142560 

15 5642.31 

-37 4915.5 

TT 

150 ±30 

4000 

3.7+^-9 

17t,p 

44 

45 

HD 143006 
HD 144432 

15 58 36.91 
160657.95 

-22 5715.2 
-27 43 09.7 

TT 

HAe 

145 ± 15 
160 ±29 

5884 

7500 

o 4-t-0.8 

nf 

70+1 

18,19,20,q 
4,e 

46 

HD 144668 

1608 34.28 

-390618.3 

HAe 

142 ± 27 

8200 

15,e 

47 

V2246 Oph 

162603.02 

-2423 36.0 

TT 

120 ±5 

5248 


21 ,r 

48 

HBC639 

162623.35 

-242059.7 

TT 

120 ±5 

5250 


21 ,s 

49 

Elias 2-24 

162624.07 

-241613.5 

TT 

120 ±5 

4266 

1-84) 

21 ,t 

50 

Elias 2-28 

162658.44 

-2445 31.8 

TT 

120 ±5 

4169 

0 

21 ,t 

51 

Elias 2-30 

1627 10.27 

-24 19 12.7 

TT/HAe 

120 ±5 

5950 

6 

5^ 

110 +'*'^ 

21 ,s 

52 

V2129 Oph 

1627 40.27 

-2422 04.1 

TT 

120 ±5 

3981 

21 ,t 

53 

V2062 Oph 

163133.46 

-2427 37.2 

TT 

120 ±5 

4900 

21 ,u 

54 

HD 150193 

164017.92 

-23 5345.1 

HAe 

216 ±68 

9500 

4,e 

55 

AKSco 

165444.84 

-3653 18.5 

HAe 

103 ±21 

6500 

4,v 

56 

KKOph 

17 1008.13 

-27 15 18.8 

HAe 

160 ± 30 

8030 

u.z. 29 

5,22,n 

57 

58 

51 Oph 

HD 163296 

17 3124.95 
175621.28 

-23 5745.5 
-215721.8 

HAe 

HAe 

124 ±4 

119 ± 11 

10000 

9200 

240)“ 

TQ + 1"6 

4,h 

4,d 

59 

HD 169142 

18 2429.77 

-294649.3 

HAe 

145 ± 15 

7500 

11+4 

21000 +^'™ 

8300+™ 

20 ,d 

60 

MWC 297 

18 27 39.52 

-03 49 52.0 

HBe 

250 ± 50 

25400 

23t,w 

61 

MWC 300 

18 2925.69 

-0604 37.2 

HBe 

1800 + 200 

19000 

24,x 

62 

RCrA 

190153.65 

-3657 07.8 

HAe 

130 ± 20 

11100 


25,26,y 

63 

64 

TCrA 

HD 179218 

190158.77 
1911 11.25 

-365749.9 
+ 154715.6 

HAe 

HAe 

130 ± 20 
254 ± 38 

6900 

9640 

4 7+^+' 
100 +^^ 

25,26,y 

4,d 


References. Distance references (t 


(5)|Hillenbrand et aLH1992); (6) 

Bally et al.H2000^; (7)|Alecian et al.H2013); (8) Kraus et al.H2013|; (9) Damiani et al.H1994|; (10) Close et al.|i 1997^; (11) Millour| 

et al. P009); (12)|Whittet et al. 

(19971; (13)|Verhoeff et al. (20101; (14)|Torres et al. (2008); (15)|Miiller et al. (20111; (16)|Preibisch & Mamajek| (20081; (17) 

ComeronlUOOSl; (18)|de Zeeuw et al. 1(19991; (19)|Bouwman et al.|(20081; (20)|Honda et al.|(2012|; (21)|Loinard et al.|(2008|; (22)|Leinert et al.|(2004); (23)|Drew| 

et al. 1^1997); (24)|Miroshnichenko et al.|^2004); (25) 

Marraco & Rydgren|(19811; (26)|Prato et al.|(20031. Effective-temperature references; (a)|Brown et al.|(20081; 

(b)|Hemandez et al.H2004^; (c)|Kenyon & Hartmann 

41995J; (d)|Eolsom et al.H2012); (e)|Alecian et al.H2013(; (f)|Mora et al.|(12002); (g)|Gray et al.H2006); (h)|van| 

den Ancker et al. ^1998); (i)|Testi et al. ^1998); (j)|Hillenbrand et al. jl992); (k)|Eolsom et al.H2008J; (l)|KrausU2009); (m) LuhmanH2007(; (n) Acke et al.H2010); 


I Prato et at. 

|et al.lpOOf T (v)' |Alenc'ar et ari^2003^ ; (w) |Acke e't al.H2008) ; (x) lMiro^'niche'nko et al.H2004^ ; (y) |Bibo et ar]jl992^ . 


A third class of objects is the TTauri stars. Full sample pa¬ 
pers about MIDI observations for these objects are still under¬ 
way. To extend our range in probed luminosities, we included 22 
representatives of the TTauri class. The sources belong to clas¬ 
sical low-mass star formation regions in Ophiuchus, Taurus, and 
Chamaeleon. 

In Table [T] an overview of the sample is given, including 
relevant stellar properties. An overview of the MIDI data that are 
used in this work is given in Appendix [B| As can be seen from 
this table, the number of observations per target varies from 1 
to 34 with a median of 6. In total, data obtained from about 240 
nights are included in the paper. The data were obtained between 
2003 and 2014, which is the full operational period of MIDI. 
A substantial amount of data are guaranteed time observations. 
Observations were done in different modes on both the 8.2-m 
Unit Telescopes (UTs) and 1.8-m Auxiliary Telescopes (ATs). 


2.3. Data reduction 


The large amount of observational data was reduced using the 
2.0 version of the EWS software package (Jaffe 2004| l, released in 
October 2012. EWS is based on a coherent integration of the ob¬ 
served visibility signal. The data reduction consists of two prin¬ 
cipal steps: (1), the extraction of the raw data from the observa¬ 
tions and (2) the calibration of these data. Step 1 is straightfor¬ 
ward and is done observation by observation (for science targets 
and calibrators). Step 2, the calibration, requires the combination 
of science and calibrator data so is more involved. 
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We analyzed the MIDI data in the form of correlated fluxes 
Fcon, which we compared to the total fluxes F of the disks 
(i.e., the spectra). Using correlated spectra is equivalent to using 
visibilities, which are the classical interferometric observables, 
but this has the advantage that the calibration does not invoke 
the photometric observations. The latter observations are often 
found to have low quality: 


1. The atmospheric and instrumental background contribution 
in the mid-infrared is strong and variable. Unlike correlated 
flux measurements, for which fringe scanning provides di¬ 
rect background subtraction at a high frequency (the frame 
rate, up to 160 Hz), photometric observations require chop¬ 
ping, which is done at a much lower frequency (2 Hz). 

2. The light path due to chopping differs for target and sky 
frames, and accordingly the background subtraction is more 
difficult. For the fringe measurement itself, the light path re¬ 
mains identical throughout the scans. 

In addition, for AT observations, the primary beam almost Alls 
the entire instrumental held of view, which complicates the esti¬ 
mation of the sky signal next to the source, on the detector. This 
can make AT photometry unusable, even for sources as bright as 
50 Jy. 

The calibrated correlated flux is calculated from the raw 
correlated flux Ccon, as follows: 


F 


corr,v 


r 


where 




pc yc 

V T V 


( 1 ) 
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(The index c denotes calibrator quantities.) Here, is the cali¬ 
brator visibility, as calculated from its known apparent diameter, 
and F‘^ is the known calibrator spectrum. The total flux F of a 
science target is calibrated in the same way as in Eq. 0 with 

= 1. Each observation therefore leads to two observables: a 
(baseline-dependent) correlated flux measurement Fcon and a to¬ 
tal flux measurement F. Einally, for a collection of observations 
of the same target, one can average the total flux measurements 
to a single high-quality spectrum. Indeed, the total spectrum is 
the equivalent of the correlated flux at zero baseline. 

The variable atmospheric transmission and coherence losses 
make the calibration in Eq. Q non-trivial: called the trans¬ 

fer function, is a time- and airmass-dependent quantity. We cor¬ 
rect for the time dependency by doing a linear interpolation in 
time of the transfer function, calculated for different calibra¬ 
tors observed the same night. The airmass dependency is cor¬ 
rected for by fitting a line to the log vs. airmass diagram 
for the entire data set and applying the corresponding airmass- 
correction factor to the observations (see [van Boekel et al.|2005| 
and|Burtscher et al.|2012|l. The uncertainty on includes both 
the intrinsic uncertainty of (photon noise) and the inter¬ 
calibrator variations of (airmass and temporal variations, 
evaluated as the standard deviation of all transfer-function ob¬ 
servations). 


All data presented in this paper will be made available in the 
Optical interferometry DataBase (OiDB), managed by the Jean- 
Marie Mariotti Center (JMMC)0 


3. Model geometry 


The large difference in number of observations (i.e., UV points) 
per target makes the model choice difficult. Generally, two op¬ 
tions can be considered for modeling the disk geometry from 
interferometric data: (1), a geometric model with a prescribed 
intensity distribution for the disk emission, and (2), a radiative- 
transfer model with a self-consistently calculated intensity dis¬ 
tribution. Given the large sample, computational reasons make 
the second approach cumbersome, and we restrict ourselves to a 
geometric model. 

In the near-infrared, inclined ring models have been success¬ 
fully used for fitting the interferometric data of disks and deter¬ 
mining their orientation (e.g., [Eisner et al.|2003 |2004[ |Benisty| 


et al. 201 1[ Kreplin et al. 2013| l. The situation in the mid-infrared 
is different. Unlike in the near-infrared, where essentially all disk 
emission comes from a compact region around the dust sublima¬ 
tion radius, the mid-infrared intensity distribution corresponds 
to a relatively extended region at different temperatures. Sim¬ 
ple geometric models (e.g., rings, uniform disks, Gaussians) do 
not represent this complexity properly. A different but related 
problem is that the estimation of disk orientations tends to be 
difficult. Examples can be found in the literature: [Kraus et al.| 
( |2008| l indicate that their disk orientation derived from near- and 
mid-infrared interferometry disagree; [Fedele et al.| ( |2008] l find a 
disk orientation for HD 135344 B that differs substantially from 
more precise estimates ( Lyo et al.|201 11 and their orientation for 
HD 101412 (i = 80 + 7°) is rather unlikely for an unobscured 
Herbig star. 

We propose the following, semi-physical model for the mid- 
infrared emission of the disks. The mid-infrared intensity dis¬ 
tribution is assumed to come from an (vertically) optically thin 
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surface layer of the disk: 


Idi^UR) = r, ByiT{R)), 


where 


T(R) = T’.ub 



( 2 ) 


In this equation, Tjub and /?sub are the sublimation temperature 
and radius of the dust, physically corresponding to the inner rim 
of the disk. The radial range of the model is R^ub < R < Rout- 
The model has two free parameters: a constant optical depth Ty 
and the temperature gradient q. Other parameters are fixed or 
precalculated: Tsub = 1500 K, /?out = 300 au, and 


J^sub ~ R-k 


/Teff 


r T* 1 

Wsub 




(3) 


(see, e.g., jPullemond & Monnier|[2OT0t note that this expres¬ 
sion assumes that the rim itself is optically thick in the radial 
direction), where J?*, Teif, and L* are the stellar radius, effective 
temperature, and luminosity, respectively. 

We assume a pole-on orientation for all disks, similar to other 
interferometric surveys based on a small amount of data per tar¬ 
get (e.g., [Monnier & Millan-Gabet|2002HMonnier et al.|2005| l. 
On the one hand, this assumption allows us to model sources 
with only a few observations (i.e., where the UV coverage in¬ 
trinsically does not allow determining the orientation). On the 
other hand, we prefer taking the same orientation for all objects 
rather than fitting orientations based on the data themselves]^ 
Since most disks have mid-infrared emission features and not 
absorption features, an orientation close to edge-on is also very 
unlikely. Under the approximation of a pole-on orientation, the 
correlated flux of the disk measured at the wavelength A — cjv 
on a projected baseline length B becomes 


Ecorr. disk,v(B) — y 


dR R ByiT(R)) Jo{2nR B/{d T)) 

C2^dRRByiTiR)) 


(4) 


where d is the distance, and Jq the 0-th order Bessel function. 
The free parameter Ty in Eq. Q is now absorbed in the free pa¬ 
rameter Edisk.v, the total disk flux at the given wavelength. The 
correlated flux F^.o„ y for an observation is then simply (the ab¬ 
solute value of) the sum of Econ. disk,v (the correlated flux of the 
disk) and ^ (the unresolved stellar flux). 

For all objects, our model is fit to the data at 4 = 10.7 /rm, 
the reference wavelength of other mid-infrared high-angular- 
resolution disk surveys ( Monnier et al.|2009{ Boley et al.|20T3] l. 
The algorithm of jForeman-Mackey et al. ( 2013| l is used for the 
parameter estimation. 


4. Results 

In Fig. the resulting fits of the temperature-gradient model to 
the data are shown per target. The corresponding parameters of 

^ As already mentioned, constraining disk orientations from mid- 
infrared interferometry turns out to be difficult. Tbe main reason for 
tbis is tbe typically low and/or non-unifonn (often unidirectional) UV 
sampling of tbe interferometric observations, in combination with an 
analytic description of tbe disk’s brightness. Interestingly, when sharp 
radial edges are present in the resolution range of MIDI, strong con¬ 
straints on the disk orientation are imprinted in the visibilities, providing 
the necessary diagnostics to fit this orientation (e.g., [Matter et al.|2014| 
[Hillen et al.|2015] Menu et al., in preparation). For continuous disks, 
a good azimuthal sampling is a necessary prerequisite. In Sect. |6.4| it 
is shown that our approximation of a pole-on orientation has no major 
influence on our results. 
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the models are shown in Table In general, the temperature- 
gradient model in Eq. Q leads to a good qualitative and quan¬ 
titative reproduction of the observed trends. The majority of the 
objects indicate a smoothly declining correlated-flux profile, a 
behavior that is captured well by the model. For a minority of 
objects, the resulting fits are rather poor, and we discuss possible 
limitations of the model below. 

First, we can consider the approximation that disks are ori¬ 
ented pole-on. Several objects with multiple observations along 
different baseline orientations indicate that the disk orientation 
indeed seems not too far from pole-on. Examples are HD 98922, 
VY Mon, and HD 142527. For other objects, such as HD 50138, 
RY Tau, and UX Ori, we see clear hints of different geometric ex¬ 
tents along different baseline angles. A likely explanation then is 
that we see these disks under a significant inclination. For the 
first two of these objects, the inclination is indeed observed to be 
relatively high (i = 56 + 4° and 66 ± 2°, see [Borges Fernandes 
et al.|20ir||Isella et al. 20101, and the inclination is also expected 
to be high for UXOri ( Grinin et al.|1994) l. Still, it is comforting 
to see that objects like HD 142666 and HD 144432, which are 
physically similar but probably have significantly different incli¬ 
nations ( Bogaert & Waelkens||1991) , are found to have similar 
extents. This suggests that the assumption of a pole-on orienta¬ 
tion has no major influence on the outcome (discussed further in 
Sect. [6^ . 

Other fits of objects indicate intrinsic discrepancies from 
the proposed disk geometry. Notable examples are 51 Oph and 
MWC 297. A closer inspection of the fits shows that the depar¬ 
ture is at least partly related to how close the inner rim can get to 
the central star. In our disk model, it is assumed that no emission 
comes from R < Rsub, where R^ub is calculated following Eq. 0 - 
The assumption puts a limit on how compact the disk emission 
can be, which translates into an upper limit on the correlated- 
flux profile. It is clear that the two mentioned objects are fit by 
a model with too many visibility nulls, which physically cor¬ 
responds to model disks that have an inner radius that is too 
extended to reproduce the compact emission. Interestingly, this 
suggests that the considered objects have a strong mid-infrared 
contribution from within the dust sublimation radius (see, e.g., 
|Acke et al.|2008| ). 

Finally, there are three objects with a mid-infrared contri¬ 
bution that is almost entirely stellar: ySPic, HD 109085, and 
V2246 0ph. As is clear from the model plots, all but the 0- 
baseline fluxes are consistent with the predicted photospherical 
fluxes. In other words, any mid-infrared excess is only detected 
in the total flux observations, and the excess is resolved out in 
the interferometric observations. The first two objects are well- 
known debris-disk objects (e.g, [Smith & Terrile |1984| |Wyatt| 
let al.||20051l, and the MIDI data give lower bounds on the dust 


location (see Smith et al.|2009 for HD 109085). V22460ph, on 
the other hand, may even lack a disk, as pointed out by jJensenj 
jet al.j ( [2009 j l. The fits of these three objects are poorly con¬ 
strained. 


5. Mid-infrared size-luminosity relation 

The temperature-gradient model in Sect.|^is expressed in terms 
of two free parameters: the total disk flux Fdisk.y and the temper¬ 
ature gradient q. The physical relevance of the first parameter is 
clear, but the second parameter needs some more attention. 

A typical temperature profile T cx is found in classi¬ 

cal models for the structure of flaring protoplanetary disks (e.g., 
Kenyon & Hartmann||1987 1. A different, limiting geometry for 
disks is a geometrically thin, passive irradiated disk, where a 
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temperature profile T cc can be derived, see, for example, 
Armitage ( 2007| l. As mentioned by the latter author, this is prob¬ 
ably the steepest profile one could expect from a passive disk. 

The temperature gradients we derive span a much wider 
range. It is therefore possibly incorrect to interpret the parameter 
as a real temperature gradient within the disk. On the one hand, 
the model description is too simplistic to capture the actual inten¬ 
sity distribution of the disk. On the other hand, the temperature 
gradient is the only model parameter for (implicitly) describing 
the geometry, so it may capture much more than the dust tem¬ 
perature alone. 

An alternative interpretation would be to see q acting as a 
radial brightness parameter. Under this point of view, we define 
a half-light radius hlr, as follows: 


U disk,v 


'>hlr 


r 

I dR2;rR/disk,v(^)- 

'-'^sub 


(5) 


In words, hlr is the disk radius within which half of the mid- 
infrared flux is emitted. Since disks may have gaps, which result 
in a discontinuity in the radial brightness profile, the use of a 
half-light radius has a clearer interpretation than a parameter that 
imposes continuity. 

The interest in a size parameter, like hlr, is that we can com¬ 
pare it directly with stellar parameters. In near-infrared interfer¬ 
ometry, it is common to compare disk sizes with the stellar lu¬ 
minosity (see jMillan-Gabet et al.|2007] and references therein). 
The size-luminosity relation we derive from our mid-infrared in¬ 
terferometric data is shown in Fig.[^ Across a luminosity range 
of almost seven orders of magnitude, we see a gradual increase 
in the size of the targets, with considerable scatter]^ 

The majority of the objects lie within the range defined by the 
expected radii at which gray, optically thin dust between 900 K 
and 250 K would be located. A handful of objects seems to be 
“oversized” with respect to this broad range. Three peculiar ob¬ 
jects in the plot are HD 109085 (#37), V22460ph (#47), and 
P Pic (#16), which have very poor constraints on their size. As 
pointed out in Sect. these objects only show a mid-infrared 
excess in the total flux, and the exact “location” of their disks 
is poorly constrained: based on the data, only a lower limit on 
the half-light radius can be derivedj^A priori, it could also have 
been expected that these debris-disk (or maybe even diskless, 
for V2246 Oph) objects are atypical. Interestingly, several other 
oversized objects are well known (pre-)transitional disks (e.g.. 


HD 169142 (#59), |Honda et al.|2012| HD 10054 6 (#34), |Benist^ 
jet al.||20T0} GM Aur (#6), Hughes et al.|2009] l, which supports 
the actual relevance of our size parameter. 

Monnier et al. ( 2009|l studied a sample of, on average, more 


massive protoplanetary disk objects using mid-infrared sparse- 
aperture interferometry on a 10 m class telescope. Although a 
different model geometry was used (a ring model), the represen¬ 
tative size range they report is very similar (our choice to show 
the 900-K and 250-K lines follows their work). Since typical vis¬ 
ibilities of the single-telescope observations are above 0.7 (or, 
> 0.5), the specific model plays a minor role, and a direct 
comparison of the size ranges is justified. We thus conclude that, 
on both large scales (probed by the sparse-aperture interferome¬ 
try) and small scales (probed by the MIDI data), the sources look 
similar. 


We should note that the term “size” only points to an apparent size: 
intrinsically, the models all have the same outer radius (300 au, see 
Sect.|^. 

^ The upper limit on hlr is put by the assumed model, which has an 
outer radius of 300 au. 
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Fig. 2. Fits of the temperature-gradient model to the 10.7 yum correlated fluxes. The gray lines represent a range of possible models for the data, 
and the dashed line is the unresolved stellar contribution at 10.7 The half-light radius Mr (in mas) corresponding to the models is included in 
the plot. 
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Fig. 2. continued. 


As already pointed out by |Monnier et al.| ( |20091 l, the rather 
poor correlation between stellar luminosity and mid-infrared 
disk size strongly contrasts with the tight correlation found 
between luminosity and near-infrared disk sizes ([Monnier & 


Millan-Gabetl |200^ [Monnier et al.|[^05[ [Millan-Gabet et ak 


2007j l. In the near-infrared, the radiation is dominated by dust 
near the sublimation temperature, located at the inner rim of the 
disk. For a given sublimation temperature, the location of this 
sublimation radius is (to first order) entirely determined by the 


stellar luminosity (Eq. Q), from which the tight correlation im¬ 
mediately follows. The region responsible for the mid-infrared 
emission, on the other hand, is much more extended. The scatter 
Fig. [3 then indicates a large variety in disk geometries. [Mon-j 
nier et al.| (2009|l mention a few physical differences that may 


explain the different appearance of disks around similar central 
stars: radial variations in grain distribution, grain growth and set¬ 
tling, and binarity. We already mentioned above that several of 
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Table 2. Fit parameters ^ and q for the temperature gradient model, per object, and the derived half-light radii hlr (in mas and an; the latter 
value involves the distance uncertainty) at d = 10.7 yum. Error bars denote the 68-% confidence intervals. 


Tdisk,F 

(Jy) 


Ur 

(mas) 


hlr 

(au) 


Tdi.sk,V 

(Jy) 


, +0.03 

0 

u.HU-t 0.003 

> 2.0 

1 5+°-3 

-0 ® 
U.J I 0 

**■■*'-0.01 
1 io+o°6 

u.j^^_0.002 

0 57+°-°* 
**■■*'- 0.01 
0 71+0.01 
**■'*- 0.01 

**■**'- 0.02 

*^■■^■^-0.02 

** ** - 0.01 

■ 3.0 0.2 

-0 6+° ° 
**■**- 1.0 

> 1.4 

**■****- 0.02 
n 44 + 0.02 
**■^^- 0.01 
0 42+0-02 

05d8® 

0 67+°-^^ 

^•^/_Qn3 

0 67+0 01 
-0.01 

^•^•^-o.os 
0 54+0 01 
-0.01 
j 9+0.1 

0 96+^^^'* 
-Q.Q3 

0 41+0.02 
**■+*- 0.01 

> 3.0 

0.568!°:™ 

0.84!°;°* 

0 49+0-02 
055^8® 

[\ f\i 


hlr 

(mas) 


Q 7+o.y 

^■^-0.8 
193 9 + 9.4 
■*■*^- 8.9 

<0.9 
6 0+°^^ 
(I2:5t,) 

14 4+0-9 

o 7+0.2 
-0.2 
19 1+^-'^ 

18.4!+* 

o 7+0.3 
-0.2 
A 7+0.4 
^••^-0.4 

9 . 3 +*^* 


hlr 

(au) 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 


LkHo- 330 
V892Tau 
RYTau 
LkCa 15 
DRTau 
GM Aur 
AB Aur 
SUAur 
HD 31648 
UXOri 
HD 36112 
HD 36917 
CQTau 
V1247 0ri 
HD 38120 
/JPic 
HD 45677 
VYMon 
HD 259431 
RMon 
HD 50138 
HD 72106 
HD 87643 
CRCha 
HD 95881 
DI Cha 
HD 97048 
HP Cha 
FMCha 
WWCha 
CV Cha 
HD 98922 


0 44+«-°2 

**■++-0,02 

28.8!°f 

1 ^+0.02 

0 73+^-^-3 
^■'■^-0.03 
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10.6+“ 

0 

21 . 2 +“ 

1 . 5 +“ 

33 

34 

35 
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TCha 

0 29+^-^^ 

^■^^-0.07 

1 7 c+0.05 

0 45k» 
^■^'^-0.02 

115+^:^ 

1 

^•^- 0.6 

36 

HD 104237 

> 1.6 

< 0.7 

< 0.10 

37 

HD 109085 

0 6+S-5’" 

^■'^'^-0.01 

90.4!3[° 

12.8!4^’ 

38 

DKCha 

19 3+^^ 

0.67!°:°! 

19.3!“ 

7 7 + 0.4 
^•'-0.4 

39 

HD 135344 B 

3 5+®-^ 
■^•*^-0.4 

0 . 68 !°;“ 

7 4 + 0.7 
'•^- 0.7 

^•^-0.4 

40 

HD 139614 

g 4 + 0.8 
^•^- 0.8 

> 1.4 

< 3.6 

< 0.5 

41 

HD 142666 


-lO.I 


1.7! 

6.4+ 


^■^-0.1 

5 9+°^^ 

-0.4 

ri lo+O.fe 
^■^°-0.02 
g 7+0.5 

0 62+^^^5 
^■^"^-0.04 

138.8+1^ 
26 9+^'^ 

^^•^-0.4 

6 30^^ ^ 

U.DU-Q , 
20 1+^^ 
574+?:^4 
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"^•^-0.2 

226.3!°3 
1 l+°^‘' 
5.! 


1 . 8 ; 

3.5; 

8 . 6 ; 

2 7 

^•/-O 

s+l 

24.8+ 


: 8 :? 

:8:i 

:83 


0 61+°°* 
**■***- 0.01 

0 60+°°* 
U.DU 0 01 
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Q 1 +0.4 
^•^-0.4 
12 

< 3.3 


A 7+0.5 

34783 

■*■+-0.9 

< 1.3 


0.48!“! 

n 77+0-06 
^■'■^-0.05 

28.9!“ 

7 7+0.7 
■^•■^-0.7 

7 7+0.1* 
■^•■^-0.6 

1 3+®-^ 
^•■^- 0.2 

0 47+0-01 
^■^'- 0.01 

33.9!^;° 

12.7!|| 

-0 6+^ ‘ 

0 71+0-01 

1 77 + 0.02 
^■^■^-0.02 

1 48+o°8 

o:8+fi! 

(11.8A) 

25.6!l;g 

229.7+1’!’, 

7 1+S 
'•^-1.9 

8 

°-^-o.i 

6 7+2-5 
-2-5 

4 2 +°-^ 

-OJ 

13.0!“ 

7 0 +O .4 

0.97^ 

19.4!ii 


0 50+^^^ 
**■■***- 0.01 

20.2![j 

5 6++“ 
■*■**-1.8 

0 72+0 0* 
u.'4 0.01 

37.3!|-3 

56.0+*** 

0.83!°:*] 

0 61+° °^ 
******- 0.02 

c 7+0.6 
^•^-0.6 

> 1.3 

< 3.4 

< 0.6 


42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 


HD 142527 
HD 142560 
HD 143006 
HD 144432 
HD 144668 
V2246 0ph 
HBC639 
Elias 2-24 
Elias 2-28 
Elias 2-30 
V2129 0ph 
V2062 Oph 
HD 150193 
AKSco 
KKOph 
51 Oph 


0.55+“( 

14.1+1* 

2 3 +“ 

58 

HD 163296 

0 56+8:81 

0 52+° °* 
^■“^^-0.01 

26.5!i| 
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59 

60 

61 

HD 169142 
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0.67+“? 
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62 
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63 
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^•■^-4 7. 

64 
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^^•^-0.4 
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1 
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-OA 

20.1 +35 
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-1.5 
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-0.5 
4 cc+0.05 
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- 0.7 

9.7!3‘ 

9 1++'' 
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12 l+°^° 
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^^•^-0.8 
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-0.3 
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27.7+H 
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■ -O.L 
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< 0.10 

0 69+0 0'7 
-0.07 
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^•^- 0.3 
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26T83 
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2 
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Notes. Numbers between parentheses need to be multiplied by 10-’ 


Table 3. Classification of the Herbig Ae stars in the sample, following 
the|Meeus et al.|p001^ classification. 


group sources 

la V892 Tau, AB Aur, HD 36112, CQ Tau, HD 38120, 
HD 259431, HD 100546, HD 139614, HD 142527, 
RCrA,TCrA, HD 179218 


Ib V1247 0ri, HD97048, HD 100453, HD 135344B, 
Elias 2-30, HD 169142 


Ila RYTau, SUAur, HD 31648, UXOri, HD 36917W , 

HD 72106, HD 95881, HD 98922, HD 104237, 

HD 142666, HD 144432, HD 144668, HD 150193, 
AKSco, KKOph, 51 Oph, HD 163296 


HD 36917 may lack any silicate feature and be a group Ilb 
source, according to Juhasz et al. (20101. Our MIDI spectrum 
for this target points to a 'weak silicate feature. 


the apparently over-sized objects are transition disks. In Sect.j^ 
we investigate the origin of the different disk sizes. 

In conclusion, our mid-infrared size-luminosity relation ex¬ 
tends the existing size-luminosity relations by an order of mag¬ 
nitude down in minimal size, toward the low-mass star formation 
domain (for the massive star regime, see Grellmann et al.|2011 
and |Boley et al.|2013| l. 


6. Size-color relations for Herbig Ae disks 

Pre-main sequence stars with similar luminosities turn out to 
have disks that are different. If all objects have the same age, 
the difference is a manifestation of the intrinsic diversity of the 
disks (e.g., in mass, viscosity, angular momentum). Conversely, 
if the objects have different ages, the diversity may hint at dif¬ 
ferent evolutionary stages. The diversity in disk sizes probably 
points to a combination of the two and forms the topic of this 
section. 


6.1. Herbig Ae disks 


More than half of the objects in our sample are categorized as 
Herbig Ae stars, and within our sample this may be the only 
group of objects with enough representatives for relevant com¬ 
parisons. We limit our further analysis to this group of objects. 
The interest in this group arises from the quoted evolutionary 
scenarios presented in Sect. and Fig. T]^ 


The Herbig Ae stars can be classified following the Meeus 
|et al. ( |2001 1 scheme: I = the “flaring” disks, and II = the “flat” 
disks. An index “a” and “b” is used to indicate whether a 10-yum 


° Although the original works propose their scenarios as representative 
for Herbig Ae/Be stars, their actual source lists are mostly limited to 
Herbig Ae objects. 
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luminosity (L^) 


Fig. 3. MIDI size-luminosity diagram for protoplanetary disks: plot of the half-light radii of the disks at 10.7 pm vs. the stellar luminosity of the 
sample targets. The gray dashed lines indicate the expected distance at which gray, optically thin dust at the indicated temperature would be located 
(cf. |Monnier et al.|2009^ . 


Table 4. Parameters of the radiative-transfer model population. 


parameter 

value 

star 

mass/lum./eff. temp. 

(1.7 Mo, 11 Lo,7000K), 

(M*, L*,7’eff) 

(2.0Mo, 21Lo, 8500K), 

(2.6 Mo, 60 Lo, 10000 K), 

(3.5 Mo, 200 Lo, 13000K) 

distance d 

140 pc 

disk 

surface density power p 

1.0 

dust mass Mdust 

O 

1 

W 

0 ^ 

disk inner radius 

R,ub, 1,2.5, 4, 6, 8, 10, 15, 20au 

disk outer radius Rout 

300 au 

minimal grain size Amin 

0.01,0.03,0.1,0.3, 1.0, 3.0yum 

maximal grain size Omax 

1 mm 

grain size distribution 

-3.5 

dust composition 

amorphous MgFeSiOd (90 %) 
carbonaceous grains (10%) 

dust-to-gas ratio 

0.01 

turbulence parameter a 

lo-^ io-\ 10-^ 10-2 

inclination i 

10, 25,40, 60, 70° 

halo 

halo inner radius 

^sub 

halo outer radius 

l-3/?sub 

dust composition halo 

same as disk, a = 0.3 - l.O/tm 

halo opacity Ty 

0.2 

interferometric model observations 

baseline lengths 

10, 20, 30,..., 130, 140m 

baseline angles 

0, 15, 30,..., 150, 165° 


6.2. Size-color diagram 

In Fig. we compare the half-light radius to the MIDI 8 - 
13 color, calculated as -2.5 log(Fv_ 8 ^m/^v,i 3 /rm)- This color is an 
estimate of the continuum spectral slope in the N band. To com¬ 
pare the objects, which have different luminosities, we normalize 
the half-light radii by the expected scaling between charac¬ 
teristic sizes and luminosities (Fig.[^. We note that the quantities 
in the plot are distance-independent]^ As is clear, sources with 
a blue A-band color tend to be small, while large objects have 
a red color, on average. Adding a color to the size comparison 
thus seems to partly decouple the degeneracy between sizes and 
luminosities. 

[Leinert et al.| ( |2004[ ) derive the mid-infrared sizes of seven 
Herbig Ae disks using MIDI interferometry and notice a simi¬ 
lar connection with the mid-infrared disk color. These authors 
interpret the correlation as evidence for the distinction of group 
I and group II by flaring versus flat disks: flaring disks will ap¬ 
pear larger, and the colder outer emitting disk surface will make 
the disk redder. The interpretation follows the classical idea of 
|Meeus et al.|P00l] l. 

In the size-luminosity diagram (Fig. |^, it is made clear that 
some objects are oversized. We pointed out that several of these 
sources are confirmed to have a gap (Sect. [^. The presence of 
radial gaps in the disk structure is not included in the flared- 
versus flat-disk picture and can be expected to alter the size-color 
relation. 

6.3. Radiative-transfer models 

To disentangle the multiple effects that play a role in our size- 
color diagram, we simulated a population of disk models with 
various parameters. We used the Monte Carlo radiative-transfer 
code MCMax ( |Min et al.|2009[ ) for this purpose. The code solves 
the temperature structure and vertical hydrostatic equilibrium 
in an azimuthally symmetric disk. Additionally, a grain-settling 
mechanism based on the strength of turbulence is implemented. 


silicate feature is seen or not. The classification for the Herbig ’ estimates for hlr and scale linear with the distance, and 

Ae stars in the sample is shown in Table colors are obviously independent of distance. 
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MIDI 8-13 


Fig. 4. Size-color diagram for the Herbig Ae stars within the sam¬ 
ple: plot of the half-light radius {Mr) of the disks, normalized to the 
luminosity versus the A-band continuum slope (calculated as 
-2.5 lo g(Fv, 8 )jm/F’v,i 3 )jm)- T he colors of the points and errorbars indi¬ 
cate the|Meeus et al.|(|2001[> classification of the sources. 


parametrized by the o'-parameter (for classical a disks). Lower¬ 
ing a decouples the larger grains and settles them to the mid¬ 
plane. 

The assumed model is a protoplanetary disk with a vary¬ 
ing inner radius and a hxed outer radius (Aout = 300 au, i.e., 
like the temperature gradient in Sect. [^, consisting of amor¬ 
phous silicates of olivine stoichiometry and carbonaceous grains 
( [Dorschner et al. 1995[ Preibisch et al.||l993| respectively). 
Of the potentially high number of parameters that can be var¬ 
ied, we limited ourselves to varying the parameters that have 
a strong influence on simulating the A-band interferometry of 
disks. The other parameters were hxed to representative values 
for Herbig Ae stars (e.g., [Mulders & Dominik|2012| l. Standard 
radiative-transfer models typically have problems reproducing 


the strong near-infrared excess of Herbig Ae stars. Mulders & 


[Dominik ( |2012| l discuss different methods for adapting standard 
disk models in order to reproduce this near-infrared bump and 
propose the use of a compact halo component in addition to the 
disk. As these authors note, the compact halo is a parametriza- 
tion of emission from the inner regions, including dust or gas 
within the inner rim. Although the near-infrared excess is not 
modeled here, the contribution of this component extends to the 
N band. We include a similar halo component here, whose prop¬ 
erties are hxed. For models without a gap, the inner part of the 
disk overlaps with the halo, and for models with a gap, the gap 
is represented by the (dust-free) region in between the halo and 
the disk. 

An overview of the parameter grid is given in Table|^ For the 
model observations, we calculated a densely sampled UV plane 
(i.e., combination of baseline lengths and angles). 


6.4. Model half-light radii: inclination and baseline effects 

To place the radiative-transfer disk models into Fig. and com¬ 
pare them to the data, we simulated MIDI observations of the 
radiative-transfer models and then ht the same temperature- 
gradient disk model (Sect. as we applied it to the data. The 
model interferometry thus allows us to derive half-light radii 
for the radiative-transfer models. Assuming that the parameter 


space for the radiative-transfer models is representative for the 
observed objects, two additional aspects in obtaining valid model 
comparisons are important: disk inclinations and baseline con- 
hgurations. 


Inclination. The temperature-gradient model does not include 
inclination effects, i.e., the model is applied under the assump¬ 
tion that the observed disks (Sect. and radiative-transfer mod¬ 
els (this section) have a pole-on orientation. In reality, protoplan¬ 
etary disks obviously have some inclination (though the Herbig 
Ae stars in our sample will not have a high inclination, since they 
are selected to have a low optical extinction). However, because 
the assumption of a pole-on orientation is applied to both the 
simulated observations of the radiative-transfer models and to 
the actual observations of the real disks, this does not introduce 
a bias in the comparison of the half-light radii. Furthermore, the 
half-light radius determined using a pole-on approximation is a 
robust parameter, and does not significantly depend on the in¬ 
trinsic disk orientation (see Appendix |A|). 


Baseline configuration. A related problem is the choice of the 
“model baseline conhguration”. The radial intensity prohle of 
the temperature-gradient model will generally be different from 
that of the actual disk. This leads to half-light radii that depend 
on the specihc baseline configuration (UV coverage). For exam¬ 
ple, the size estimate of an inclined disk along the major and 
minor axes differs. For each radiative-transfer disk model there 
is therefore a corresponding range of half-light radii. 

To calculate this range of half-light radii for a given 
radiative-transfer model, we fit the temperature-gradient model 
to each individual UV point. Given that the full UV coverage 
of a single model contains 168 points (14 baseline lengths x 12 
baseline angles), we thus get the same number of half-light radii, 
per model. The plots in the next section will show the range of 
the half-light radii we get. 

The set of MIDI observations defines a baseline configura¬ 
tion for each target. Except for four objects, all Herbig Ae stars 
have at least two UV points. Having a UV coverage with several 
baseline lengths and baseline angles will lead to a more con¬ 
strained value for the half-light radius, which averages over the 
size estimates of observations on a single baseline. Even if the 
temperature-gradient description is only approximative and as¬ 
sumes a pole-on orientation, our calculation of the model half- 
light radii, a parameter that is robust to inclination effects, repre¬ 
sents a very conservative basis for comparing with the observed 
half-light radii. 

6.5. Theoretical versus observed size-color relations 

Eigure|^shows the size-color diagram for the model population. 
The hgure is split up into panels that show the four different stel¬ 
lar luminosities for which we have made radiative-transfer mod¬ 
els. Eor all models, the range of half-light radii (see Sect. |6.4| i is 
indicated. The axes ranges are identical to those in Eig.|^ 

Eor all luminosity cases, the colors and sizes of models with¬ 
out a gap (i.e., with Ri^ = Rsub) spans a similar range (a broad¬ 
ening downward is seen for the highest luminosity case). This 
range is relatively compact: gapless disks can only have a lim¬ 
ited range of half-light radii and mid-infrared colors. Introducing 
a gap in the disk by moving out the inner radius of the disk com¬ 
ponent substantially expands the possible range of colors and 
half-light radii. This occurs in a qualitatively similar fashion for 
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Fig. 5. Size-color diagram for the population of radiative-transfer models for the four different luminosity cases (11 Lq, 21 Lq, 60 Lq, and 200 Lq). 
The vertical lines indicate the range of sizes that can be obtained for the individual models. The 5 % smallest and 5 % largest mid-infrared sizes 
are shown in transparent colors, giving an idea of the typical (90 %) range of model sizes. In black, the models without a gap are shown. The other 
colors show models with a gap between the halo and the disk components (for an increasing disk inner radius 


all luminosities. As can be expected, large gaps will move the 
half-light radius outward, and the larger emitting surface of rel¬ 
atively low-temperature material makes the disk redder. Similar 
gap sizes correspond to bluer disk colors for more luminous cen¬ 
tral objects. This is a simple consequence of the disk wall being 
warmer for the more luminous central objects. 

The model size-color plot allows us now to interpret the 
MIDI size-color diagram for Herbig Ae stars. In Fig.|^ we com¬ 
bined the observational diagram (Fig. with the approximate 
model ranges in Fig. Except for a few outliers, which are all 
in group Ib, the range defined by the observations largely cor¬ 
responds to the model range. Focusing first on the models with 
a continuous disk (gray zone in Fig. we notice that only a 
limited number of sample stars have the size-color combination 
falling into this area of the diagram. The majority of these targets 
are group II sources. Of the group la group, most objects have 
characteristics corresponding to the gapped sources. They share 
this property with a significant number of the group II sources. 

[Maaskant et al.|p0I3| l suggest all group I sources may have 
gaps, whereas it is proposed that group II sources have contin¬ 
uous, gapless disks. Our size-color diagram not only confirms 
that most group I objects have (inner) disks that differ from stan¬ 
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dard, continuous disks. It also shows in a simple way that there 
may be a significant overlap in structure of the inner disk region 
for both groups or, at least, for group la and group II objects. A 
comparison of the overlap region with the model diagrams even 
suggests that the structure of these group II disks differs from 
that of continuous disks, and they may have gaps. Gaps in group 
II disks would alter the picture that group I and group 11 disks 
form two strictly separated populations. The above findings are 
discussed further in Sect. |7] 


6.6. Non-thermal mid-infrared emission 


Objects that clearly are not represented by our model population 
are the group Ib disks with very blue colors. Group Ib objects 
lack the A-band silicate feature by definition, a property that is 
associated with their large gaps in the radial dust distribution 
(e.g., [Honda et al.||2012j [Maaskant et al.|[2013[ [Carmona et al.| 

2014[ l. Instead of the silicate feature, their A-band spectra typi¬ 
cally have strong emission peaks of polycyclic aromatic hydro¬ 
carbons (PAHs; Juhasz et al.|2010| . These peaks are intimately 
related to the gap nature: the decrease in dust emission makes it 
easier to observe the PAH emission ([Maaskant et al.|2013[), and 
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Fig. 6. Merged size-color diagram, combining observations (Fig. 
with the model population (Fig. |^. In gray, we indicate the range of 
models without a gap. In yellow, we show the approximate range of 
models with a gap between the inner halo and the disk component. We 
note that both regions indeed overlap. Group Ib sources with colors de¬ 
ferring from the model population are V1247 Ori (#14), HD 135344 B 
(#39), Elias 2-30 (#51), and HD 169142 (#59). 


the ionization imprint in the PAH emission traces the density 
throughout the disk ( Maaskant et al.|2014| l. 

The emission mechanism of PAHs is based on the electronic 
excitation by ultraviolet photons and the subsequent emission in 
molecular stretching and bending modes. In (and just outside) 
the N band, emission peaks of these modes are centered at 7.7, 
8.6, 11.3, and 12.7 /rm. The transient heating mechanism, a prop¬ 
erty that PAHs share with very small grains (VSGs), is a poten¬ 
tial source of mid-infrared emission coming from the outer disk 
regions. 

The effect of PAHs and VSGs on A-band interferometry is 
twofold. First, a strong PAH contribution, notably around the 
typically strong and wide 8//m band, makes the MIDI 8-13 
color bluer. Second, a strong mid-infrared contribution of PAHs 
or VSGs in the outer disk will make a disk appear much larger 
than if its emission were thermal. The relative contribution of 
these molecules or grains to the N band is the strongest for group 
Ib objects, which explains why they have deviant disk colors and 
sizes in Fig.[^ In a sense, the connection between apparent disk 
size and color, usually coupled by the grain temperature of parti¬ 
cles in thermal equilibrium with the radiation field, is decoupled 
by the stochastically heated grains. 


6.7. Limitations of the modei popuiation 

A direct comparison of the diagram suggests that the model pop¬ 
ulation covers the necessary range of parameters to interpret our 
observations. However, we need to be careful with the interpre¬ 
tation. 

The model population consists of two-component models 
with a halo and a disk, which can be decoupled by a gap. This 
model is representative for some Herbig Ae disks, but certainly 
not for all. Notable exceptions are the transitional disks that 
still have a massive inner disk. Example objects are AB Aur (la, 
[Honda et al.||M0l l and HD 97048 (Ib, [Maaskant et ar||20T3l l. 
Including these more complex disk structures is, at least in prin¬ 
ciple, possible in the disk population, but introduces several new 


parameters. The latter would increase the number of models, 
hence the model computation time, significantly. The extra mod¬ 
els also will not give any new insight; the observational diagram 
is already fully covered by the models (apart from some group 
Ib objects), which indicates that new models will only point to 
model degeneracies. 

This point brings us to the main limitation of our compari¬ 
son; mid-infrared sizes and colors are (obviously) insufficient for 
assessing the full disk structure. The following problems can be 
identified; 


A fixed dust composition; our models use a fixed dust com¬ 
position. Varying the silicate composition and/or the contin¬ 
uum opacities may slightly change the MIDI 8-13 color. 

A fixed halo component; the halo component is kept fixed for 
the model population. Making the halo component weaker 
will make the near-infrared excess disappear, which is ex¬ 
cluded by typical Herbig Ae SEDs ([Mulders & Dominik| 


20121. Increasing this component makes the mid-infrared 


emission bluer and more compact. 

Gaps in the outer disk; the mid-infrared wavelength range is 
only sufficient for exploring gaps in the inner disks. Gaps in 
the outer disk will only be visible at longer wavelengths. In 
the absence of any other gaps, a disk will appear as standard 
in the N band. 

Gap sizes; the gap size in a given part of the model diagram 
depends on the luminosity of the central star. When translat¬ 
ing the color and size of an object into a possible gap size, 
one needs to take the luminosity of the central star into ac¬ 
count. 

Model size range; the range of sizes (or half-light radii) that 
can correspond to a given color seems much larger from the 
gapped-disk models than what is actually observed. This is 
mainly a model bias. As explained in Sect. 6.4 each disk 


model has a maximum range of model sizes. Having a few 
observations, which is typically the case for our MIDI data, 
strongly confines the size estimate of a disk (Sect. [6^. 


In summary, our model population might partly underestimate 
the range in colors that disks can have, and it also has the ten¬ 
dency to overestimate the typical size range of disks. 

These limitations of the model population may hamper a di¬ 
rect comparison with the MIDI size-color diagram. However, 
it is clear that most effects are not changing our main conclu¬ 
sions in Sects. [63] and [6.6[ some group II disks are incompatible 
with standard disks, and the contribution of stochastically heated 
grains is needed to interpret the diagram. 


7. Tracing the evolution of Herbig Ae disks 

Ideally, understanding the structure of Herbig Ae disks in the 
mid-infrared demands a detailed individual modeling of a large 
sample of targets with a high number of observations at high spa¬ 
tial resolution. A large number of objects have indeed been ob¬ 
served with mid-infrared interferometry, but the number of ob¬ 
servations per target is generally too low for constraining sophis¬ 
ticated models. We modeled the existing, archival observations 
in a simple way and interpret the results in this section. 

7.1. Overiap in spatiai properties group I and II 

Protoplanetary disks in hydrostatic equilibrium that are irradi¬ 
ated by a central star will have a limited range of possible mid- 
infrared colors and (apparent) sizes, determined by the repro¬ 
cessing of the stellar radiation by the dust. Under the reasonable 
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assumption that the objects are largely azimuthally symmetric, 
any deviation from the standard color-size relation must point to 
differences in the radial and/or vertical structure of the disk. 

The idea that most group I objects may possess large ra¬ 
dial gaps is almost exclusively based on observational evidence 
from the outer disks around Herbig Ae stars ( Honda et al.|2012] 
Maaskant et al.|[2013|l. Our mid-infrared interferometric analy¬ 


sis of disks provides more insight into the nature of the inner 
disks of group I objects. From our comparison between the ob¬ 
served sizes and colors of Herbig Ae disks and sizes and col¬ 
ors of radiative-transfer models in Sect. we found that most 
group I objects cannot be interpreted as having “standard” con¬ 
tinuous disks. They fall in a region where the radiative-transfer 
models have gaps in their radial dust distribution. This direct ev¬ 
idence for gaps in the inner disks of group I objects endorses the 
idea that most of these disks are transitional. Whether all group 
I disks in our sample have gaps cannot be confirmed, since the 
observations of at least a few objects (e.g., RCra (#62)) are in 
agreement with there being a continuous inner disk. 

For the group II objects. Fig. indicated that several have 
colors and sizes that agree with models for continuous disks. 
This agrees with the classical picture that group II disks are 
gapless and flat. However, a number of group II sources differ 
from “standard” continuous disks. They fall in an overlapping re¬ 
gion with the gapped group I sources. Also the radiative-transfer 
models in the overlap regions possess gaps in their radial dust 
distribution. Is it correct to interpret these group II sources as 
having a gap? For at least two group II disk sources, gaps have 
already been claimed. Based on near- and mid-infrared interfer¬ 
ometric observations (IOTA and MIDI, respectively) |Schegerer| 
et al. (|2013 i modeled the disk around HD 142666 as having a gap 


with a radius of 0.8 au, plus an additional inner disk between 0.3 
and 0.35 au. |Chen et al.| ( [20T^ , on the other hand, observed and 
analyzed the disk around HD 144432, and their interferometric 
observations (AMBER and MIDI) point to a similar model with 
a gap radius of 1 au and a compact 0.2-au radius inner ring. The 
two sources are group II sources in the overlap region between 
group I and II sources (#41 and #45), and both sources fall out¬ 
side of the range of disk models without gaps. 

In addition to the above examples, millimeter observations 
of the group II source RY Tau indicate that the inner disk envi¬ 
ronment may be dust-depleted ( |Isella et al.||2010| l. The lack of 
millimeter emission for this object may not immediately imply a 
physically gapped disk (e.g., |Dong et al.|2012 1 . Still, in the size- 
color diagram in Fig. RY Tau (#3) lies at the top of the region 
that is compatible with disk models without gaps. This object 
may therefore be an example of a gapped group II disk that is 
missed by the simple size and color criteria. Our conservative 
way of estimating the possible range of half-light radii may in 
fact impede the identification of RY Tau as a gapped source, and 
even more group II objects with large half-light radii may have 
gaps. 

We conclude that, although the direct translation of sizes and 
colors in the presence of gaps in group II disks might not be 
straightforward, there seems to be a promising correspondence 
between the two. The discovery that also some group II disks 
have gaps alters the existing picture of group I and II disks. 


7.2. The observational picture of Herbig Ae disks 

With MIDI we observe in the 8-13 micron atmospheric window 
{N band) and probe mainly warm dusty material at T > 200 K. 
This limits the radial range that we see to several tens of au from 
the star, because at larger radii the dust is not warm enough to 
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/MIDI range x MIDI range 



Fig. 7. Comparison between the spatial range probed by MIDI (in terms 
of its spatial resolution and the temperature domain) and the possible, 
observed disk structures. We distinguish between group I and group 
II disks, with or without silicate feature (a vs. b). The green and red 
colors indicate whether MIDI provides the spatial information to fully 
characterize the presence and onset (inward) of possible disk gaps. 


emit significantly in this spectral region]^ The range of spatial 
scales that is sampled by our observations and in which we can 
directly detect the signature of a disk gap is 1 to 15 au in the best 
cases, but depends on the UV coverage. 

Figure |7] gives an overview of possible disk geometries for 
Herbig Ae stars. We distinguish between the following cases: 

1. Groups I and 11 ( |Meeus et al.||20()T] ). These indicate the 
amount of far-infrared excess emission and are thought to 
physically correspond to objects with flared (group I) or flat 
(group II) outer-disk geometries. 

2. Subgroups a and b ( |Meeus et al.||200l] l, indicating whether 
emission features from silicate dust are seen (a) or not (b). 
Physically this corresponds to a configuration where small 
(< few micron) silicate grains are present at the disk sur¬ 
face in the right temperature range to emit in the N band 
(subgroup a), or not (subgroup b). Following the [Maaskant 
|et al.| ( |2013| l picture, group Ib objects have very large gaps 
in the region which otherwise dominates the silicate emis¬ 
sion. Transitional group la objects still have strong silicate 
emission, which is dominated by either the outer disk (la in 
Fig.[7|ii) or the inner disk (la in Fig.jT^i). 


® The radial range seen in the 10-pm spectral region depends on the 
stellar luminosity (scaling with These considerations only apply 
to dust that is in thermal equilibrium with the radiation field, such as 
typical silicate dust. In contrast, particles that are stochastically heated, 
such as PAHs, can be observed to much greater distances of up to sev¬ 
eral 100 au in some HAe stars. Such particles are small enough to reach 
temperatures of several 100 K upon absorption of a single optical/UV 
photon, making them “light up” briefly at mid-infrared wavelengths. 
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The above classification does not provide a sufficient characteri¬ 
zation for fully interpreting our MIDI results. In addition to this 
commonly used framework, we include a third dimension; 


3. The spatial range of a possible disk gap. There may either 
be no gap, a gap within the spatial range probed by MIDI, 
or a gap outside the spatial range probed by MIDI. For the 
three geometries on the left in Fig. [7] MIDI has the capacity 
to correctly infer whether a gap is present in the disk and to 
measure the (inward) onset of the gap(s). For the right three 
geometries, either the main gap is too extended to be identi¬ 
fied with MIDI (alone), or it lies outside the temperature or 
resolution range that is probed. In the group II sources, with 
a disk gap, the spatial extent of the main gap is smaller than 
in group I sources. 


The geometries depicted in Fig.[7]correspond to models that have 
been proposed for Herbig Ae stars. Two exceptions can be iden¬ 
tified: the group Ila disk with a gap far from the star and the 
group la disk with two gaps (Figs. |7]/ an d|7]/i, respectively). 

The first case, the group Ila source with a gapped outer disk 
may currently be missing owing to an observational bias; the gap 
in such objects is difficult to detect because the radial width is too 
small for the gap to have a significant effect on the overall SED. 
In addition, group II disks are faint in scattered light (e.g., Garufi 
|et al.|2014| l, and the gap lies outside the disk region probed with 
MIDI and can only be detected with long-baseline sub-mm/mm 
observations with ALMA or the VLA. 

The second case, the geometry for a double-gapped group 
la disk in Fig. |7] requires some additional discussion. Several 
group la disks are characterized by large gaps of tens of au in 
size. Examples are AB Aur, HD 36112, and HD 142527 with in¬ 
ner radii of the outer disks of 88 au, 73 au, and 130 au, respec¬ 
tively ( [Honda et al.||2010{ [Andrews et al.||2OTT| [Verhoeff et al.[ 
2011[l. The physical extent of these gaps, and the fact that they 


still show strong 10-/tm silicate emission (i.e., subgroup a), in¬ 
dicates that their 10-//m emission is coming exclusively from an 
inner disk, a component that is potentially small and partially 
depleted of dust (e.g., [Andrews et al.|2011[ l. Mid-infrared inter¬ 
ferometry of such a source will point to a compact disk with 
small sizes and blue colors similar to a strongly settled contin¬ 
uous disk. In other words, the sources are expected in the range 
of models without gaps in Eig. Instead, sources like AB Aur 
(#7), HD 36112 (#11), and HD 142527 (#42) have sizes and col¬ 
ors that agree with models with small inner gaps. A possible in¬ 
terpretation is that these group la sources have additional small 
gaps in the inner disk regions, similar to the newly identified 
gapped group II sources. The depicted geometry in Eig. |7] cor¬ 
responds to such a geometry. Alternatively, the inner disk ge¬ 
ometry might be altered in a different way and contain an ex¬ 
tended halo component, as was proposed for HD 142527 ( [Ver-[ 
[hoeff et al.|2011[ we note that such a component strongly differs 
from the compact tenuous halo of our radiative transfer models 
in Sect. [6.3| l. Recent work by [Marino et al.[ ( |2015[ ) suggests that 
the inner disk of HD 142527 is actually strongly inclined with 
respect to the outer disk. 


7 . 3 . Effects of gaps on other disk observables 

Is the formation of a gap at a specific location in the disk the 
fundamental characteristic that determines the observational ap¬ 
pearance of Herbig Ae disks? Eorming a gap by, for instance, 
the formation of a sufficiently massive planet will lower the lo¬ 
cal surface density and largely isolate the inner disk from the rest 



Fig. 8. Near- and mid-infrared excess SED of the group Ila stars in 
the sample. The SEDs are normalized to the stellar flux in the J band 
(T = 1.24pm). The full yellow lines indicate sources with observed 
or suspected gaps. The shaded green zone show the 25 % to 75 % per¬ 
centile of the other objects. Sources with gaps generally have slightly 
lower near-infrared excesses than gapless sources. 


of the disk. If this process occurs close to the star, continuous ac¬ 
cretion of the inner-disk material can partly drain out the inner 
disk. A lower amount of hot material will lead to a decrease in 
near-infrared radiation, an effect that will be visible in the ob¬ 
ject’s SED. 

Eigure gives an overview of the SEDs of the group II ob¬ 
jects in the sample, where the stellar photosphere has been sub¬ 
tracted. In order to compare objects of different luminosities and 
at different distances, we normalized the fluxes to the photo- 
spheric contribution in the J band (d = 1.24pm). The aim of 
the figure is to compare the following classes of sources: 

- group Ila objects with observed/expected gaps within 15 au 
(based on Eig. |^, i.e., sources HD 72106 (#22), KKOph 
(#56), HD 142666 (#41), UXOri (#10), HD 144432 (#45), 
and AKSco (#55); 

- other group Ila objects. 


Objects with an inner gap tend to have slightly lower near- 
infrared fluxes than the other group members P The mid-infrared 
fluxes, on the other hand, are similar to the other targets within 
the group. The lower near-infrared excesses for sources with in¬ 
ner gaps is likely related to the partial accretion (“drain-out”) of 
the inner disk, as alluded to above. 


Acke et al. (20091 find a strong anti-correlation between the 


7-pm excess and the [30/13.5] continuum flux ratio for Her¬ 
big Ae/Be stars. In their interpretation, a higher inner disk re¬ 
sults in a stronger shadowing of the outer disk, hence a bluer 
mid-infrared disk color. Conversely, the shadowing effect of a 
lower inner-disk rim is lower, so these disks appear redder. The 
relatively low near-infrared excess for disks with an inner gap 
now provides an alternative explanation for this correlation. In 


® A notable exception is the upper group Ila source with a (suspected) 
inner gap. This source is KKOph (#56). |&eplin et aL] ( [2013^ model its 
near- and mid-infrared interferometric data with a large envelope with 
polar cavities, in addition to a (continuous) disk. This geometry deviates 
from a standard disk, and possibly explains why the source is atypical 
in our sample. 
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an evolutionary scenario driven by gap formation, objects with a 
strong near-infrared excess and blue colors have the most primi¬ 
tive inner disks, where a strongly optically thick inner rim casts 
a shadow on the outer disk. Objects where gap formation in the 
inner disk has led to a partial drain-out have a lower shadowing 
effect. These disks are less optically thick, have a more illumi¬ 
nated outer disk, and will appear redder. 

An important diagnostic for the disk structure is the PAH 
emission. Based on the work of |Acke et al.| ( [20T0| l, |Maaskant| 
et al. ( 2014|l compile a list of 48 Herbig Ae/Be sources and com- 


Observations that spatially resolve Herbig Ae disks chal¬ 


pare their PAH luminosity to the |Meeus et al.| ( |2001| l classifica¬ 
tion. Whereas the group Ib objects all are found to have a rel¬ 
atively strong PAH emission, the PAH contribution in group la 
and Ila objects is not very different. For a similar fraction of 
sources (6/13 for group la and 12/28 for group Ila), PAH emis¬ 
sion is not detected, and the intrinsic emission strengths overlap 
(with a wider spread for the Ila sources). The strong contribution 
for PAHs in group Ib sources can be attributed to their very large 
gaps, which lower the continuum emission and increases the ir¬ 
radiation to UV emission ( |Maaskant et al.|[2013| |2014| l. Group 
la and Ila sources have either smaller or no gaps, which largely 
excludes this effect. Similarly, the stronger molecular gas emis¬ 
sion for group I sources than for group II sources can be related 
to a more progressed dust clearing within the disks of the former 
group, resulting in a more efficient heating of the outer layers 
( |Meeus et al.|2013] l. 

7.4. Evolutionary classification for TTauri and Herbig stars 
Herbig Ae stars are the more massive counterparts of the more 


numerous TTauri stars. Spitzer surveys of clusters (e.g., Lada 
e^^lj |20()^ [Currie et al.||2009l |Merm et al.||2010] jCieza et'aL 


2012 jl have yielded new insight into the characteristics and evo¬ 


lution of TTauri disks. The SEDs of TTauri stars have been 
classified in four groups: primordial disks, disks with large in¬ 
ner holes (or gaps), homologously depleted (or anemic/weak- 
excess) disks with weak overall infrared excesses, and debris 
disks. The different classes of objects in between primordial and 
debris disks led to the idea that different evolutionary channels 


exist toward debris-disk or diskless objects (e.g., Lada et al. 
2006 jl, governed by the different dispersal processes in proto¬ 


planetary disks (e. g., jCieza et al. |2012| l. An evolutionary sce¬ 
nario proposed by [Currie et al. ( 2()09[ see also |Curn^|2010| l is 
that primordial disks evolve into either homologously depleted 
disks through simultaneous dust clearing at a wide range of disk 
locations or disks that developed large inner holes (or gaps) as 
part of an inside-out clearing process. The Maaskant et al. ( [2013| l 
scenario can be interpreted as the Herbig Ae equivalent for the 
above evolutionary scenario, where group II disks are the ho¬ 
mologously depleted disks and group I disks the disks with large 
holes or gaps. 

The classification and evolutionary scenarios of T Tauri stars 
are essentially based on SEDs, i.e., spatially unresolved data. 
Even though radiative-transfer modeling may partly lift degen¬ 
eracies based on simple color criteria (e.g., [Merln et al.|[25T0l 
[Currie & Sicilia- Aguilar|201 l[ l and may correctly help to iden¬ 
tify large inner holes, the structural information in SEDs is still 
limited. The group I and group II SED classification for Herbig 
Ae stars is an excellent example of this. Initial ideas suggested 
grain growth and settling as the main mechanism to transform 
group I into group II disks ( [Meeus et al.||200T] [DuIIemond &| 
Dominik|2004a|b| l, and several works on grain growth endorsed 


lenged this evolutionary picture. The realization of Maaskant 
et al. P013| l that all group I objects may have gaps shows that 
many “primordial” group I disks are in fact in a far-evolved stage 
of gap formation. Our results provide additional evidence for the 
gapped nature of group I disks. The presence of large gaps sug¬ 
gests that many disks bear massive, or multiple (e.g., |Zhu et al.j 
201 1[ [Dodson-Robinson & Salyk 201 l| l, planets. If all group 

I objects turn out to have gaps]™ it may well be that group I 
objects indeed form a separate evolutionary channel ( Maaskant] 
et al.[2013 see Sect. [^. A possible common ancestor for Herbig 
Ae disks, like the primordial disk for T Tauri stars, has not (yet) 
been identified. 

Our spatially resolved observations of a large sample of Her¬ 
big Ae disks, identifying a significant number of gapped group 

II disks, now provide a second alteration of existing evolution¬ 
ary scenarios. The possibility that these objects, in a far stage 
of dust growth and settling, may also be in an advanced stage 
of planet formation and disk dissipation is intriguing in its own 
right. Erom an evolutionary point of view, the new population of 
gapped disks may point to new links between the existing classes 
of objects. 

7.5. Evolutionary implications of gapped group II disks 

A full target-by-target characterization, including the MIDI data 
and other high-angular-resolution constraints, is an important 
next step to be taken for understanding the newly identified pop¬ 
ulation of gapped group II disks. In the current picture, the gaps 
in flat disks are small and are located in the inner few au. Elat 
disks with small gaps in the outer disk, as in Eig. [7]/, have not 
(yet) been discovered. Based on the limited but qualitatively new 
information given by the MIDI observations, we now discuss 
possible evolutionary implications of this newly identified pop¬ 
ulation. 

A population of gapped group II disks provides direct evi¬ 
dence that multiple disk-dispersal mechanisms - grain growth 
and gap formation - are contributing together to the evolution of 
individual disks. The classification of objects based on their spa¬ 
tially unresolved appearance (e.g., the decision tree in Eig. 10 
Williams & Cieza]|2011l is therefore a simplification of the 


actual dispersal process. 

Different ideas for incorporating the new class of objects 
in current evolutionary scenarios can be thought of Eirst, the 
gapped group II objects may be evolved group I objects, with a 
continued grain growth and settling of the outer disk. Contrary 
to the gaps discovered in group I disks, gaps in group II disks are 
small. This scenario would therefore require the collapse of the 
outer disk and the closing of large gaps to be coeval processes. 
Though this cannot be excluded, it is not likely to be a common 
phenomenon. In addition, group Ib objects, which are character¬ 
ized by extreme gap sizes and a lack of small silicate grains in 
the inner several tens of au, most probably do not evolve into 
quasi-continuous disks with a clear silicate feature. Two other 
evolutionary ideas are the following: 

1. Even though group II objects have a grain population that 
experienced significant evolution, gaps may still be develop¬ 
ing, yet the gaps are smaller than for the group I objects, and 
group I and II sources present a parallel evolution. 


these evolutionary ideas (e.g., Acke et al.|2004[ Bouwman et al. 


[2008i|M^s et al.|2009[[Juhasz et al.|2010[l. 


Other recent and upcoming works show additional evidence of 
gapped group I objects (e.g., [Matter et al.|2014[|van der Plas et al.[2015[ 
Khalafinejad et al., in prep.; Menu et al., in prep.). 
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Fig. 9. Possible evolutionary connections between Herbig Ae/Be sources. Three populations of objects are identified: flat continuous disks (group II 
sources), flat disks with gaps (group II sources), and flaring transition disks (group I sources). Gapped group II sources are most likely descendants 
of group II sources without gaps. The arrows indicate two possible evolutionary ideas for the ancestors of transitional group I sources: green arrow 
= influx from a common ancestor; yellow arrow = influx from group II sources with gaps (cf. Fig.[^. 


2. Group II objects may be the precursors of (some) group I 
objects. The link between the two could be giant planet for¬ 
mation, leading to the carving of a gap and the dynamical 
excitation of the disk material. 

In both scenarios, group II objects without gaps are the precur¬ 
sors of group II objects with gaps. The main difference between 
the two ideas is the origin of group I objects. As we depict in 
Fig. 1^ (green arrow), the first scenario assumes that the group 
I population represents the outcome of an isolated evolutionary 
channel. The second scenario (yellow arrow if Fig. [^, on the 
other hand, assumes at least a partial influx from group II ob¬ 
jects. 

The first scenario can be embedded naturally in a common- 
ancestor scenario (for Herbig Ae/Be stars, Maaskant et al.|20TJ 


of a few Myr have either debris disks or no disks at all (e.g., 


for T Tauri stars, |Currie et al.|2009| ). The origin of the different 
evolutionary outcomes, in terms of the initial conditions, is still 
not understood. A common ancestor for Herbig Ae/Be stars has 
not yet been identified. It is therefore not clear whether low- and 
intermediate-mass objects follow exactly the same evolution. For 
instance, unlike for evolved T Tauri stars, a large number of the 
intrinsically evolved Herbig Ae objects are accreting (e.g., |Gar-| 
|cia Lopez et al.|[2006| l, and clear inner holes (as opposed to 
gaps) seem not to be common. Finally, Herbig Ae disks are often 
found to be relatively massive (e.g., |Acke et al.|20()4| l. Identify¬ 
ing group II disks with the population of (intrinsically low-mass) 
homologously depleted TTauri disks is therefore problematic. 

The second scenario, which incorporates a possible evolu¬ 
tionary connection between group II and group I objects, is dis¬ 
cussed in the next section. 


7.6. Can (some) flat disks evolve in flaring transition disks? 

Age information is very important for verifying different evolu¬ 
tionary ideas. Young stellar clusters, generally having well estab¬ 
lished agesj^contain young intermediate-mass objects. It turns 
out that a large number of these objects in clusters with ages 


Sicilia-Aguilar et al.|2006[ Currie et al.|2009[ Currie & Kenyon 

2009[ [Hernandez et al.|2009| ). It is clear that Ae Herbig Ae stars 
with massive disks may therefore be a fairly selective population 
of young stars. For about half of the objects in our sample, an 
age estimate can be found in [van Boekel et ar| ( |2005| ). The age 
values show that many objects in our sample are relatively old 
(> 2 Myr), and no clear evidence for a different age distribution 
for group I and II objects can be inferred. These age estimates of 
isolated objects have large uncertainties, however. Arguing for 
different scenarios based on the individual age estimates for the 
objects is therefore intricate, and the reasoning below is based 
on more heuristic arguments. 

The initial grain growth in disks can be very fast. Within on 
the order of 0.1 Myr, the entire submicron grain population of a 
disk will experience growth and grains with sizes up to 100 ;um 
or 1 cm can be formed (e.g., |Birnstiel et al.|2010 and references 
therein). Naturally, by the time the sources have cleared most 
of the envelope and become optically visible (Lada) class II ob¬ 
jects, their disks may have a settled grain population and tend 
to be “flat” (i.e., Meeus group II without gap). Depending on 
the overall grain growth, disks will be less or more settled (e.g., 
jjuhasz et al.|2010) l. 

The growth from cm-sized pebbles to km-sized planetesi- 
mals and finally planets is a more complex process, involving 
the overcoming of several growth barriers (for a recent review, 
see jJohansen et al.|2014| l. The typical timescale for the formation 
of a massive planet, which is needed for carving out a clear gap 
in the disk, remains an important unknown quantity. Depending 
on the exact grain growth models (compact versus fluffy aggre¬ 
gation), the expected timescales for the transition from cm-siz ed 
grains to 100-m sized bodies i s 10^ yr ( Okuzumi et al.|2012| to 
10® yr ( Windmark et al.|[2012 i, at a radial distance of 1 au (see 
also Schmitt et al.|1997 1 . Growth to genuine planetesimals pos¬ 
sibly even requires mechanisms that lead to a large-scale con- 


'* However, as stated by Bell et al. |2013[ , the ages of star-forming 
regions may be significantly underestimated. We refer to|Cloutier et al.| 


( |2014| ) for a further discussion of the cluster populations in terms of the 
previous and new cluster ages. 
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centration of particles (e.g., |Johansen et al.|2007| l, and the char¬ 
acterization of these mechanisms is an ongoing process. 

On the time scale of 1 - lOMyr, objects of > 10 Me can 
form, which are massive enough to keep hold of the H and He 
gas (e.g., [Pollack et al.|1996| l. These objects quickly evolve into 
young gas giant planets of > 1 Mjup, which start clearing a gap 
in the disk, but this gap is not yet very large and does not signif¬ 
icantly affect the SED. Depending on at what radius in the disk 
the planet forms and opens the gap, we may see it with MIDI if 
it occurs in the right range (i.e., group II sources with a gap that 
is detected; Fig. [^i) or may not see it if the gap lies outside of 
the range probed by MIDI (i.e., group II sources with a gap that 
is present but remains undetected; Fig.|^). 

For a sufficiently massive planet, the gap will become so 
large (group I) that a significant range of t emp eratures is “miss¬ 
ing”, yielding the typical “dip” in the SEDp^The massive body 
dynamically excites the population of planetesimals (e.g., |Krijt| 
& Dominik 201 1[), leading to a collisional cascade and a renewed 


production of (sub-)micron sized grains (e.g., [Bouwman et al. 
|2003| l. This increases the opacity of the disk, hence its ability 
to flare ( [Dullemond & Dominik|2004a| l. Furthermore, the direct 
dynamical interaction of the disk with the massive planet may in¬ 
crease its scale height at the inner edge of the outer disk: energy 
dissipated in spiral shock waves heat the disk midplane (|Harker 


|& Deschi||2002| l, increasing the pressure scale height. Depend¬ 
ing on the spatial extent of the gap, we may or may not see the 
mid-infrared silicate emission. If there is still much material in 
the right temperature range, we are seeing a group la source. In 
some of these disks, the emission is dominated by material in 
the inner disk (e.g., HD 142527), in others by the material in the 
outer disk (e.g., HD 100546; see Fig. 0. If the gap becomes so 
large that the entire region responsible for the mid-infrared sil¬ 
icate emission is empty, we are seeing a group Ib source (e.g., 
HD 169142). 

The above scenario explains how group II objects may 
evolve into group I objects via group II objects with gaps. It re¬ 
mains to be shown from a modeling point of view whether such 
a transformation is possible. It may not be likely or required that 
all flaring transitional disks are descendants of flat continuous 
disks. Still, this scenario does not require any unidentified com¬ 
mon ancestor, and it naturally explains the occurrence of all the 
kinds of disks we see, including the absence of group Ilb sources. 
Indeed, for these sources to exist, very large gaps in group II 
disks would be required (for suppressing the silicate emission), 
hence very massive or multiple planets. The latter would turn the 
disk into a group I object long before the gap is large enough to 
see no silicates. 


8. Summary and conclusions 

Linking signatures for the structure of protoplanetary disks 
to their evolutionary status is one of the primary interests of 
disk observations. High-angular-resolution observations provide 
valuable constraints on this discussion. 

Mid-infrared interferometric observations of 64 protoplane¬ 
tary disks, obtained with the MIDI instrument on the VLTI, were 
collectively analyzed in this work. The typically low number of 
observations per target imposes the use of simple disk models. 
The temperature-gradient model we used leads to size estimates 
for the disk region emitting in the mid-infrared. 


Multiple planets may be required to form large gaps (e.g., Zhu et al. 
2011 Dodson-Robinson & Salyk 201 1[(. 
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For the total sample, going from low-mass T Tauri objects to 
massive YSOs, we see a gradual increase in characteristic sizes, 
but with a considerable spread. Several outliers in the trends are 
well known for having large disk gaps, implying that the sizes 
derived using simple disk models provide a diagnostic for the 
presence of disk gaps. 

To quantitatively interpret observed trends, a population of 
radiative-transfer models was assembled. Limiting the compari¬ 
son to the Herbig Ae disks, which with 35 representatives form 
the largest group of objects within the sample, leads to four im¬ 
portant findings: 


Most group I objects have sizes and colors that are compat¬ 
ible with having gapped inner disks. This endorses the idea 
of [Honda et al. ( |2012[ l and Maaskant et al. ( 2013j l that many, 
and maybe all, group I objects have gaps. 


2. Some group II disks, thought to have flat, gapless disks, are 
indeed compatible with such a configuration. 

3. Several other group II objects have sizes and colors that are 
incompatible with such a configuration. Instead, they share 
characteristics with group I objects, which are all thought to 
have gaps. They also are compatible with radiative-transfer 
models of transition disks. For some of the mentioned group 
II disks, it has already been claimed that they have gaps. 

4. Some disks have colors that differ strongly from those of 
the radiative-transfer models. All these objects are group Ib 
sources, whose mid-infrared emission (by definition) lacks 
a silicate feature. We argue that the emission of PAHs and 
VSGs alters the mid-infrared appearance of these disks and 
decouples the relation between disk color and (apparent) disk 
size, usually governed by the dust that is in thermal equilib¬ 
rium with the radiation field. 


The population of group II disks with gaps makes a strict struc¬ 
tural separation between continuous flat disks and flaring transi¬ 
tion disks obsolete. 

Gaps in massive protoplanetary disks are often associated 
with the presence of massive planets. It is remarkable that the 
new population of group II disks also shows gaps. This poten¬ 
tially indicates that these objects are also in an advanced stage 
of gas-giant planet formation. 

From an evolutionary point of view, flat disks with gaps 
(gapped group II) most likely are descendants of continuous flat 
disks (group II without gaps). This idea can naturally be embed¬ 
ded in a common-ancestor scenario where flaring transition disks 
(group I) form a separate evolutionary channel. Gap formation 
by massive planets may also have a stronger structural impact 
on the initially settled disk. Gaps growing in size may lead to a 
group I appearance. A dynamical excitation of large bodies may 
reintroduce a micron-sized grain population. In addition, planet- 
disk interaction may dynamically increase the height attained by 
dust particles. The result of these interactions is a disk that may 
appear to be group I, with large gaps and a flaring outer disk. 

Full grain-growth models for disks, which include the forma¬ 
tion of and the interaction with a planetary body, can shed new 
light on the link between observed disk structures and evolution¬ 
ary stage. 
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Appendix A: Temperature-gradient modeis: the 
poie-on approximation 

The assumption of the temperature-gradient model (Sect. that 
disks are oriented pole-on is incorrect, at least for some disks. In 
this appendix, we comment further on this approximation and its 
possible effect on the resulting half-light radii. 

The pole-on approximation is equivalent to assuming that the 
inclination i is zero, an orientation for which the disk’s posi¬ 
tion angle P.A. is not defined. For the interferometric observa¬ 
tion, this orientation has the advantage that the model is inde¬ 
pendent of the baseline angle. Indeed, the relative angle between 
the disk’s position angle and the baseline angle is what generally 
plays the role in defining the model orientation. For this reason, 
the pole-on approximation has generally been used for interfer¬ 
ometric surveys with few observations per target (e.g., |M onnier| 
|& Millan-Gabet|2QQ2[[Monnier et al.|2QQ5l t. 

Intrinsically, the pole-on approximation is only justified for 
pole-on or mildly inclined (e.g., i < 20°) disks. However, a sig¬ 
nificant number of the disks in our sample will have a stronger 
inclination. To justify the use of a pole-on disk geometry for de¬ 
termining the half-light radius of these disks, we perform the 
following s imul ation. We take two of the radiative transfer mod¬ 
els of Sect. 6.3 with the same stellar/disk parameter^but with 
two different inclinations: i - 10° (nearly pole-on) and i - 60° 
(strongly inclined, close to the maximum for a non-obscured 
central object). For each of the two models, we calculated the 
half-light radius with the pole-on temperature-gradient model, 
for a random set of five interferometric observations (i.e., five 
UV points). This experiment was repeated 500 times, and his 


tograms of the determined half-light radii are shown in Fig. A. 1 


First, the Monte Carlo simulation shows that even for this strong 
inclination difference, the median half-light radius for both dis¬ 
tributions differs by only 10%. Second, the fit of the strongly 
inclined disks is slightly biased toward underestimating the half- 
light radii found for the (almost) pole-on disk, and the range 
of possible size estimates is 20-25 % wider. These minor differ¬ 
ences allow us to conclude that the mid-infrared half-light radius 
of a pole-on temperature gradient model is a robust parameter, 
even for disks that are strongly inclined. The conclusions based 
on Fig. 1^ (for which the vertical axis is on a logarithmic scale) 
are thus unaffected by this approximation. 

Two alternatives for this pole-on approximation can be con¬ 
sidered, for which we show below that they provide less robust 


Parameter values: M^, = 2.0 Mg 
Raub, flmin = 0.01pm, a = 10“^ 


L* — 21 Lq, T^i 


8500 K, R„ 
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Fig. A.I. Results of a Monte Carlo simulation for testing the influence 
of the pole-on approximation of the temperature-gradient models on 
inclined disks. The blue and red histograms show the distribution of 
(normalized) half-light radius estimates for a (almost) pole-on disk (i = 
10°) and a strongly inclined disk (i = 60°), respectively. The median 
size estimates differ by 10 %. 



Fig. A.2. Results for the same Monte Carlo simulation as in Fig. A. 1 
(for the radiative-transfer disk with i - 60°), but with a temperature- 
gradient model that also includes the inclination / and the position angle 
P.A. as free parameters. Clearly, neither the half-light radius, nor i and 
P.A. are well constrained. The radiative-transfer disk has cos i = 0.5 
(i.e., i = 60°) and P.A. = 0°/180°. 


or less confined results. First, it is obviously possible to extend 
the fit of the temperature-gradient model to include the disk in¬ 
clination and position angle as fit parameters. We did this exper¬ 
iment for the above radiative-transfer model disk with i - 60°. 
In the first histogram in Fig. A.2 we see that the inferred half- 
light radius is much less constrained than under the pole-on ap¬ 
proximation in Fig. |A.1| The two other histograms show the in¬ 
ferred inclination and position angle, neither of which are well 
constrained. It is clear that the originally robust size parameter 
(under the pole-on approximation) is not robust when the disk 
orientation is assumed to be free. 

A second option is to fix a non-zero inclination for the 
temperature-gradient model and determine half-light radii with 
this inclined geometric model. To avoid biases related to the un¬ 
known position angle, the applied model needs to be fit at the full 
range of position angles (P.A. = 0° to 180°/360°). The result of 
such a fit is a range of half-light radii (for the varying position 
angles) rather than a single value. Part of this size range will 
come from models that are oriented perpendicularly to the ac¬ 
tual disk orientation. For strongly inclined disks, these half-light 
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radius estimates will therefore be less precise than when a pole- 
on model is taken. The result is a less confined size estimate than 
for the pole-on approximation. The conceptually easier pole-on 
approximation, which we have shown to be robust (even when 
disks are strongly inclined), was therefore the preferred approach 
in this work. 


Appendix B: Overview of MiDi observations 
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Table B.l. Overview of observations 


name 

time 

(UTC) 

base 

Bp 

(m) 

P.A. 

(°) 

time 

(UTC) 

base 

Bp 

(m) 

P.A. 

(°) 

time 

(UTC) 

base 

Bp 

(m) 

P.A. 

(°) 

LkHu 330 

2011-12-11,04:08 

U2U3 

40.5 

54.1 

2011-12-11,04:46 

U2U3 

42.9 

51.4 

2011-12-12,05:24 

U1U4 

130.1 

56.3 

V892 Tau 

2004-12-28,05:02 

U2U3 

45.6 

46.1 

2004-12-28, 05:05 

U2U3 

45.7 

45.8 

2004-12-30, 04:08 

U3U4 

48.5 

92.0 


2004-12-30,04:11 

U3U4 

48.0 

91.7 

2009-09-04, 09:43 

U1U2 

37.2 

37.5 





RYTau 

2004-11-01,04:43 

U4U2 

78.9 

96.5 

2004-11-01,04:57 

U4U2 

81.0 

94.6 

2004-11-04, 07:51 

U3U4 

48.6 

92.0 


2004-12-28,01:47 

U2U3 

32.1 

54.6 









LkCa 15 

2011-12-11,02:31 

U2U3 

29.9 

47.1 

2011-12-11,03:10 

U2U3 

33.9 

49.7 

2011-12-11,03:22 

U2U3 

35.0 

50.2 


2011-12-11,03:29 

U2U3 

35.8 

50.4 









DR Tau 

2004-11-01,06:47 

U4U2 

89.0 

83.3 

2005-01-01,02:40 

U3U4 

61.0 

106.2 

2005-01-01,02:50 

U3U4 

60.4 

105.6 

GM Aur 

2011-12-11,05:28 

U2U3 

41.7 

52.5 

2011-12-11,05:41 

U2U3 

42.5 

51.7 





AB Aur 

2004-12-28,03:50 

U2U3 

39.7 

54.2 

2004-12-28, 06:00 

U2U3 

45.9 

44.5 

2004-12-30, 05:20 

U3U4 

42.8 

85.4 


2005-01-02, 05:03 

U2U4 

83.6 

66.8 

2005-12-23,02:48 

U1U2 

33.9 

37.4 

2005-12-24, 03:03 

U3U1 

69.1 

47.5 


2005-12-26,02:47 

U1U4 

112.4 

75.0 

2005-12-26, 03:00 

U1U4 

115.3 

73.8 

2009-01-21,01:02 

HOGO 

29.5 

82.7 


2009-01-21,01:13 

HOGO 

30.0 

81.4 

2009-01-21,01:45 

HOGO 

31.1 

77.7 

2009-01-21,01:57 

HOGO 

31.4 

76.3 


2009-01-21,02:05 

HOGO 

31.6 

75.3 

2009-01-22, 02:16 

EOGO 

15.9 

73.6 

2009-01-22, 03:32 

EOHO 

47.4 

63.6 


2009-01-22, 03:39 

EOHO 

47.2 

62.6 









SU Aur 

2006-11-10,06:26 

U1U4 

120.1 

71.5 









HD 31648 

2007-02-05, 02:08 

U2U3 

42.8 

51.2 









UX Ori 

2004-12-28, 06:25 

U2U3 

44.5 

44.8 

2004-12-30, 07:02 

U3U4 

33.0 

136.0 

2007-10-25, 07:34 

U2U3 

44.9 

42.1 


2007-10-25,07:49 

U2U3 

45.5 

43.1 

2007-11-23,03:09 

U1U4 

97.7 

48.8 





HD 36112 

2004-12-30,06:12 

U3U4 

39.6 

89.2 

2005-01-02, 06:18 

U2U4 

76.7 

62.6 

2005-12-23, 03:49 

U1U2 

39.6 

36.1 


2005-12-24, 03:57 

U3U1 

77.5 

44.5 

2005-12-26, 03:54 

U1U4 

120.8 

69.8 

2007-01-03,02:13 

U2U3 

31.1 

51.3 


2007-01-06,01:49 

U1U3 

62.2 

38.7 









HD 36917 

2011-12-11,07:38 

U2U3 

45.3 

45.9 

2011-12-12, 07:34 

U1U4 

116.7 

62.9 





CQ Tau 

2005-01-01,03:38 

U3U4 

59.8 

103.5 









VI247 Ori 

2011-12-11,06:57 

U2U3 

46.6 

45.9 

2011-12-12, 06:53 

U1U4 

126.0 

62.7 

2011-12-12,07:01 

U1U4 

124.9 

62.5 

HD 38120 

2005-12-23,01:51 

U1U2 

49.3 

8.4 

2005-12-26,01:47 

U1U4 

101.9 

49.8 

2005-12-27, 04:27 

U2U3 

45.9 

43.4 

/? Pic 

2006-12-07,04:09 

U3U4 

58.4 

94.7 

2006-12-07, 04:43 

U3U4 

60.2 

101.0 

2006-12-07, 05:26 

U3U4 

61.7 

109.0 


2006-12-07, 06:28 

U3U4 

62.5 

120.7 

2007-01-03, 03:44 

U2U3 

43.7 

44.0 

2007-01-04,04:10 

U2U3 

42.7 

48.5 


2007-01-04,05:15 

U2U3 

39.5 

57.9 

2007-02-04, 02:48 

U2U3 

40.9 

54.3 





HD 45677 

2006-10-17,05:44 

HODO 

40.9 

41.1 

2006-11-14,04:26 

DOGO 

22.8 

49.0 

2006-11-14, 07:31 

DOGO 

31.9 

71.0 


2007-10-04,07:28 

EOGO 

12.1 

52.9 

2007-12-09, 04:06 

GlHO 

69.9 

175.4 

2007-12-11,08:04 

GlHO 

71.4 

19.7 


2008-01-12,07:31 

EOGO 

10.3 

78.7 









VY Mon 

2010-01-13,03:23 

GIDO 

67.2 

128.4 

2010-01-14, 03:45 

HODO 

63.2 

74.0 

2010-01-15,04:19 

EOGO 

16.0 

72.5 


2010-01-16,04:37 

EOHO 

47.8 

71.5 

2010-01-17,03:43 

HOGO 

31.8 

73.7 

2010-01-19,02:56 

KOGl 

75.1 

28.7 


2010-01-19,05:51 

KOGl 

89.2 

35.8 

2010-01-22, 03:58 

KOAO 

127.9 

72.2 





HD 259431 

2004-10-30, 08:57 

U4U2 

89.4 

82.1 

2004-11-01,05:22 

U4U2 

56.0 

90.2 

2004-12-29, 06:41 

U2U3 

46.5 

46.0 


2004-12-30,02:33 

U3U4 

59.6 

113.7 

2004-12-31,04:26 

U3U4 

61.6 

107.6 

2005-01-01,05:43 

U3U4 

54.7 

105.9 


2007-02-08, 05:41 

U1U3 

102.0 

36.8 

2007-03-10, 03:12 

U1U2 

56.4 

35.1 





R Mon 

2009-01-19,04:02 

HODO 

64.0 

72.9 

2009-01-19,04:27 

HODO 

63.8 

71.9 

2009-01-19,05:11 

HODO 

62.0 

69.6 


2009-01-21,03:15 

HOGO 

31.4 

74.0 

2009-01-21,03:26 

HOGO 

31.7 

73.8 

2009-01-22, 03:59 

EOHO 

48.0 

72.5 


2009-01-22, 04:03 

EOHO 

48.0 

72.4 

2009-01-22, 04:28 

EOHO 

47.6 

71.3 

2010-01-13, 04:29 

GIDO 

61.9 

130.9 


2010-01-14,05:05 

HODO 

63.3 

71.1 

2010-01-16, 02:24 

EOGO 

14.1 

75.1 

2010-01-16, 05:19 

EOHO 

46.7 

69.9 


2010-01-19,04:54 

KOGl 

86.1 

35.0 

2010-01-19,04:58 

KOGl 

86.4 

35.1 

2010-01-21,02:40 

KOAO 

120.6 

74.8 


2010-01-21,06:02 

AOGl 

56.8 

116.9 









HD 50138 

2007-12-09, 03:19 

GlHO 

68.9 

167.8 

2007-12-09, 03:40 

GlHO 

68.5 

169.8 

2007-12-10, 06:33 

GlHO 

68.4 

9.5 


2007-12-12,05:09 

GIDO 

71.5 

130.0 

2007-12-12, 06:00 

GlHO 

68.1 

6.9 

2007-12-13,03:10 

GIDO 

66.6 

129.9 


2007-12-26,07:45 

GIDO 

58.1 

151.3 

2008-11-09, 05:48 

EOGO 

12.3 

60.6 

2008-11-10, 05:43 

EOHO 

36.7 

60.4 


2008-11-10,05:47 

EOHO 

37.1 

61.0 

2008-12-27,06:35 

EOGO 

15.5 

74.0 

2008-12-28,01:39 

EOGO 

10.0 

50.8 


2008-12-28,02:36 

EOHO 

36.9 

60.8 

2008-12-30,01:35 

HOGO 

20.3 

51.7 

2009-01-21,06:01 

HOGO 

27.6 

74.1 


2009-01-21,06:08 

HOGO 

27.1 

74.0 

2009-03-08,01:22 

EOHO 

47.6 

73.4 





HD 72106 

2005-12-30,07:35 

U1U4 

124.2 

69.5 

2005-12-31,04:28 

U1U4 

127.4 

40.8 

2006-03-12,00:53 

U2U4 

88.2 

72.0 
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Table B.l. continued. 


name 

time 

(UTC) 

base 

Bp 

(m) 

P.A. 

(°) 

time 

(UTC) 

base 

Bp 

(m) 

P.A. 

(°) 

time 

(UTC) 

base 

Bp 

(m) 

P.A. 

D 

HD 87643 

2006-02-26, 04:33 

DOGO 

31.3 

76.4 

2006-02-28, 06:22 

AOGO 

57.1 

101.3 

2006-03-01,01:23 

AOGO 

63.5 

38.2 


2006-03-01,05:03 

AOGO 

61.1 

85.1 

2006-04-19, 04:34 

DOGO 

25.4 

123.1 

2006-05-23, 02:08 

AOGO 

51.6 

119.9 


2006-05-25,01:23 

AOGO 

54.4 

110.6 









CR Cha 

2006-05-19, 04:23 

U1U3 

57.3 

108.4 

2006-05-19, 04:43 

U1U3 

55.7 

114.0 

2007-05-06, 23:38 

U3U4 

58.7 

87.5 


2007-05-06, 23:49 

U3U4 

58.9 

90.1 

2007-06-26, 00:24 

U3U4 

61.8 

145.5 

2007-06-26, 00:40 

U3U4 

61.9 

149.4 

HD 95881 

2004-06-06,01:35 

U3U1 

66.4 

77.1 

2004-06-07,01:47 

U3U1 

64.7 

80.6 





DI Cha 

2006-05-19, 05:23 

U1U3 

53.9 

123.9 

2007-05-07, 00:33 

U3U4 

59.6 

98.7 

2007-05-07, 00:44 

U3U4 

59.8 

101.3 


2007-06-26,01:09 

U3U4 

61.9 

154.2 

2007-06-26,01:26 

U3U4 

61.9 

158.2 

2008-05-18,01:26 

U1U3 

71.2 

64.0 


2008-05-19, 00:46 

U1U3 

73.3 

56.2 









HD 97048 

2004-12-28,07:19 

U2U3 

38.9 

28.2 

2004-12-30, 07:55 

U3U4 

58.5 

84.0 

2011-03-20,04:22 

U2U3 

36.4 

59.4 


2011-03-22, 03:23 

U3U4 

59.4 

95.8 

2011-03-22, 03:27 

U3U4 

59.5 

96.8 





HP Cha 

2005-05-28, 23:52 

U1U2 

38.7 

44.8 

2005-05-30,01:01 

U3U4 

61.2 

127.2 

2005-05-30,01:13 

U3U4 

61.3 

130.1 


2008-04-24, 00:20 

U2U4 

89.3 

65.4 









FMCha 

2007-05-07,01:48 

U3U4 

60.8 

116.4 

2007-05-07,01:59 

U3U4 

60.9 

119.0 

2008-05-19,00:00 

U1U3 

76.6 

45.2 

WWCha 

2005-05-29, 00:59 

U1U2 

37.0 

57.4 

2005-05-29, 23:53 

U3U4 

60.5 

110.9 

2005-05-30, 00:05 

U3U4 

60.7 

113.8 

CV Cha 

2007-06-26, 02:24 

U3U4 

62.2 

170.9 

2007-06-26, 02:35 

U3U4 

62.2 

173.5 

2008-04-21, 04:46 

U1U3 

65.2 

83.3 

HD 98922 

2006-12-30, 08:43 

U3U4 

59.9 

99.5 

2007-01-03, 08:37 

U1U4 

128.6 

56.7 

2007-02-06, 06:02 

U1U3 

96.6 

27.1 


2008-04-08, 07:20 

EOGO 

11.1 

132.9 

2008-05-24, 00:57 

HODO 

60.9 

85.0 

2008-05-25,01:52 

HODO 

57.1 

97.1 


2008-05-25, 02:02 

HODO 

56.4 

99.2 

2008-05-25,02:14 

HODO 

55.4 

101.9 

2008-06-01,02:11 

GIDO 

71.2 

160.1 


2008-06-01,02:21 

GIDO 

71.3 

161.9 

2008-06-01,03:41 

GIDO 

71.5 

177.3 

2009-01-19,05:38 

HODO 

63.1 

42.6 


2009-01-21,04:29 

HOGO 

30.9 

28.4 

2009-01-21,05:05 

HOGO 

31.3 

36.9 

2009-01-22,05:13 

EOHO 

47.1 

39.7 


2009-04-20, 05:22 

HOGO 

25.5 

113.2 

2009-04-23, 04:09 

EOHO 

42.3 

99.3 

2009-04-23, 04:57 

EOHO 

39.2 

110.0 


2009-05-24, 23:53 

HODO 

63.2 

73.2 

2009-05-25, 00:26 

HODO 

62.1 

79.7 

2009-06-03,02:51 

EOGO 

12.3 

118.7 


2010-01-17,07:08 

EOGO 

16.0 

59.6 

2010-01-17,07:54 

EOGO 

15.9 

68.8 

2010-01-17,08:44 

EOGO 

15.6 

78.4 

HD 100453 

2004-12-28,09:01 

U2U3 

44.3 

36.8 

2005-01-02, 09:08 

U2U4 

89.4 

78.8 

2005-12-27,06:16 

U2U3 

46.1 

7.2 

HD 100546 

2003-06-17,00:11 

U3U1 

74.4 

60.5 

2004-06-03,01:24 

U3U2 

34.8 

74.2 

2004-12-28,08:08 

U2U3 

41.3 

30.8 


2005-12-27,08:06 

U2U3 

41.4 

29.5 

2006-02-13,04:01 

EOGO 

16.0 

39.2 

2006-02-16, 05:47 

EOGO 

15.8 

66.6 


2006-02-16, 08:07 

EOGO 

14.9 

99.3 









TCha 

2011-03-20, 05:13 

U2U3 

35.9 

61.0 

2011-03-23,03:23 

U2U3 

37.7 

39.1 

2011-03-24, 03:35 

U1U4 

121.3 

58.1 

HD 104237 

2005-12-27,06:53 

U2U3 

39.3 

10.5 

2006-03-12, 05:15 

U2U3 

36.9 

53.2 

2006-03-12,06:28 

U2U3 

35.2 

69.5 


2006-03-12, 06:40 

U2U3 

34.9 

72.3 

2006-05-16, 02:25 

U1U3 

71.1 

63.7 

2006-05-17,00:54 

U1U4 

119.9 

69.2 


2006-05-18, 02:09 

U2U3 

35.1 

70.6 

2010-05-04, 23:54 

KOGl 

65.6 

16.2 

2010-05-05,00:13 

KOGl 

65.3 

19.8 


2010-05-05,00:16 

KOGl 

65.2 

20.6 









HD 109085 

2007-03-06, 07:49 

U1U3 

101.9 

39.2 

2007-03-08, 05:48 

U3U4 

62.4 

109.0 

2007-03-09, 07:20 

U1U3 

102.3 

38.3 


2007-03-10, 05:27 

U2U3 

45.6 

36.8 

2007-03-10, 05:50 

U2U3 

46.1 

39.1 

2008-03-20, 07:14 

U1U2 

56.2 

34.6 


2010-01-30, 09:34 

U2U4 

87.2 

85.4 

2010-01-30, 09:42 

U2U4 

86.5 

85.9 

2010-02-03, 04:45 

U2U3 

40.8 

8.1 


2010-02-03, 04:54 

U2U3 

40.9 

9.9 

2010-02-03, 08:33 

U2U3 

46.5 

41.5 

2010-02-03, 08:41 

U2U3 

46.6 

42.2 

DK Cha 

2013-05-02, 03:02 

U1U3 

75.6 

47.9 

2013-05-03, 02:33 

U1U2 

40.1 

35.3 





HD 135344B 

2006-04-15, 05:34 

U1U4 

130.1 

53.7 

2006-05-14, 05:04 

U3U4 

62.4 

115.5 

2006-05-14, 07:39 

U3U4 

56.3 

144.6 


2006-05-16, 03:40 

U1U3 

101.8 

27.7 

2006-06-10, 23:58 

U3U4 

50.4 

87.0 

2006-06-13,23:42 

U1U2 

56.5 

4.5 


2006-06-14, 02:13 

U1U2 

55.7 

24.8 

2006-07-11,00:54 

U2U4 

89.2 

83.4 

2006-07-13,03:37 

U2U1 

47.7 

42.9 


2010-04-25, 03:58 

U2U4 

84.8 

64.3 

2010-04-25,04:11 

U2U4 

85.9 

66.5 

2010-05-28, 02:02 

U1U2 

56.4 

14.9 


2010-06-28, 00:40 

U2U4 

88.5 

73.4 

2011-03-22, 08:04 

U2U4 

89.3 

82.2 

2011-04-13,03:47 

U2U3 

46.4 

15.8 


2011-04-14, 05:21 

U1U3 

102.2 

24.3 

2011-04-16, 04:16 

U2U3 

46.5 

22.6 

2011-06-15,00:53 

U3U4 

57.8 

97.4 


2011-06-15,01:32 

U3U4 

60.5 

102.9 

2011-06-15,03:29 

U1U4 

124.0 

69.6 

2011-06-15,04:07 

U1U4 

117.9 

74.5 


2012-05-10, 06:19 

U1U3 

94.8 

43.1 

2012-06-05,23:46 

U2U4 

75.5 

46.8 

2012-06-06, 02:43 

U2U4 

89.4 

79.4 


2012-06-06, 04:32 

U2U4 

83.7 

95.9 









HD 139614 

2010-04-26, 03:16 

U3U4 

49.8 

81.7 

2010-04-26, 04:02 

U3U4 

54.8 

90.0 

2010-04-26, 05:35 

U3U4 

61.2 

105.1 


2011-04-14,07:58 

U1U2 

52.4 

34.5 

2011-04-18,03:35 

U1U2 

56.1 

2.1 

2011-04-18,05:15 

U1U2 

55.6 

16.4 

HD 142666 

2004-06-08,04:18 

U3U1 

102.2 

36.7 

2004-06-08, 06:50 

U3U1 

91.0 

43.7 

2004-06-08, 07:01 

U3U1 

89.6 

43.7 

HD 142527 

2003-06-14, 00:32 

U3U1 

102.2 

10.9 

2003-06-14, 00:54 

U3U1 

102.1 

14.3 

2005-06-23, 00:58 

U1U4 

129.6 

45.8 


2005-06-23,02:10 

U1U4 

129.7 

58.0 

2005-06-23,02:21 

U1U4 

129.3 

59.7 

2005-06-24, 06:33 

U3U4 

56.1 

159.1 


2005-06-27, 05:07 

U2U3 

36.7 

61.2 

2005-06-27,05:18 

U2U3 

35.8 

62.3 

2005-06-27,06:18 

U2U3 

30.0 

67.7 


2005-06-27, 23:54 

U2U4 

83.3 

56.4 

2005-06-28,01:09 

U2U4 

88.4 

70.8 

2005-06-28,01:20 

U2U4 

88.8 

72.6 


2012-06-06,00:16 

U2U4 

77.3 

41.9 
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Table B.l. continued. 


name 

time 

(UTC) 

base 

Bp 

(m) 

P.A. 

(°) 

time 

(UTC) 

base 

Bp 

(m) 

P.A. 

(°) 

time 

(UTC) 

base 

Bp 

(m) 

P.A. 

(°) 

HD 142560 

2005-05-26, 00:56 

U3U4 

44.2 

78.6 

2005-05-26, 03:43 

U3U4 

61.1 

104.6 

2005-05-26, 07:55 

U3U4 

55.4 

150.1 


2005-06-24, 00:11 

U3U4 

53.4 

90.4 

2005-07-23, 00:54 

U2U3 

45.5 

43.3 

2005-07-23, 04:37 

U2U3 

31.0 

63.1 


2005-08-26,02:18 

U2U4 

61.0 

125.3 

2005-08-26,01:55 

U2U4 

65.3 

119.2 

2006-05-15, 02:29 

U2U3 

46.4 

17.2 


2006-05-15,05:02 

U2U3 

46.0 

40.3 









HD 143006 

2007-06-26, 03:19 

U3U4 

61.1 

116.7 









HD 144432 

2003-06-17,02:15 

U3U1 

102.3 

26.2 

2003-06-17,04:46 

U3U1 

98.6 

41.1 

2003-06-17,04:51 

U3U1 

98.2 

41.5 


2003-06-17,05:07 

U3U1 

97.0 

42.5 

2004-02-11,08:52 

U3U2 

45.1 

18.5 

2004-02-12,09:56 

U3U2 

46.1 

29.3 


2004-04-11,05:55 

U3U2 

46.0 

28.1 

2004-04-12, 06:16 

U3U2 

46.3 

31.6 

2006-03-12,08:13 

U2U3 

46.3 

30.6 


2006-05-16, 08:41 

U1U3 

86.2 

46.1 

2006-05-16, 08:51 

U1U3 

84.6 

46.3 

2006-07-11,05:19 

U1U3 

83.2 

46.3 

HD 144668 

2004-04-10, 02:44 

U3U2 

46.3 

173.4 

2004-04-10, 05:23 

U3U2 

46.6 

21.6 

2004-04-11,05:16 

U3U2 

46.6 

21.1 


2004-04-11,09:09 

U3U2 

42.3 

52.4 

2004-04-12, 09:46 

U3U2 

39.8 

56.4 

2004-06-28, 02:06 

U3U1 

100.8 

30.9 


2004-06-28, 05:22 

U3U1 

83.3 

50.6 

2004-06-28,00:15 

U3U1 

102.4 

15.4 

2004-06-28, 04:26 

U3U1 

90.9 

46.1 


2004-09-30, 00:02 

U4U2 

63.5 

123.8 

2006-03-12, 09:39 

U2U3 

45.6 

42.1 

2006-03-12, 09:59 

U2U3 

45.1 

44.6 


2006-07-10, 23:37 

U1U4 

127.9 

41.8 

2009-06-30, 04:09 

KOAO 

114.8 

88.6 

2009-06-30, 04:20 

KOAO 

112.6 

90.3 


2009-07-01,00:55 

KOAO 

126.4 

58.9 

2009-07-01,01:42 

KOAO 

128.0 

66.7 

2009-07-01,01:44 

KOAO 

128.0 

67.2 


2009-07-01,01:47 

KOAO 

128.0 

67.6 

2009-07-01,01:50 

KOAO 

128.0 

68.0 

2009-07-01, 05:42 

AOGl 

83.4 

154.1 


2009-07-01,05:54 

AOGl 

82.8 

156.8 

2009-07-01,05:58 

AOGl 

82.7 

157.7 

2009-07-01,06:01 

AOGl 

82.5 

158.6 


2009-07-02,02:18 

KOGl 

87.9 

27.8 

2009-07-02, 02:22 

KOGl 

87.7 

28.2 

2009-07-04, 00:10 

EOHO 

46.5 

53.0 


2009-07-04,00:14 

EOHO 

46.6 

53.7 

2009-07-06, 00:30 

EOHO 

47.3 

58.1 

2009-07-06, 00:34 

EOHO 

47.4 

58.8 


2009-07-06, 03:31 

EOHO 

44.1 

86.3 

2009-07-06, 03:35 

EOHO 

43.8 

86.9 

2009-07-07, 03:53 

HOGO 

28.2 

90.3 


2009-07-07, 03:56 

HOGO 

28.0 

90.9 









V2246 Oph 

2012-05-10,07:59 

U1U3 

97.8 

42.1 

2012-06-06, 00:54 

U2U4 

67.6 

53.7 

2012-06-06,01:01 

U2U4 

68.8 

55.1 

HBC 639 

2005-04-18,04:46 

U2U4 

74.5 

61.2 

2005-04-18,04:50 

U2U4 

75.2 

61.9 

2005-04-19, 03:53 

U2U4 

65.3 

51.0 


2005-04-19,04:03 

U2U4 

67.2 

53.3 

2005-08-26, 00:20 

U3U4 

59.3 

121.6 





Elias 2-24 

2006-05-19, 08:31 

U1U3 

90.9 

44.4 

2007-05-07, 04:32 

U3U4 

58.3 

101.9 

2007-05-07, 04:39 

U3U4 

58.8 

102.5 


2007-05-07,05:19 

U3U4 

61.4 

106.0 

2007-05-07, 05:30 

U3U4 

61.8 

107.0 

2008-04-21,08:15 

U1U3 

101.3 

38.3 


2008-05-18,03:56 

U1U3 

101.1 

21.8 

2008-05-18,04:38 

U1U3 

102.0 

27.4 

2008-05-19,05:51 

U1U3 

102.2 

35.5 

Elias 2-28 

2012-06-06, 03:14 

U2U4 

87.5 

75.4 









Elias 2-30 

2012-06-06,01:30 

U2U4 

74.0 

60.7 

2012-06-06,01:41 

U2U4 

75.7 

62.5 

2012-06-06,01:47 

U2U4 

76.8 

63.6 

V2129 Oph 

2011-05-14,06:03 

U3U4 

62.3 

113.0 









V2062 Oph 

2007-06-26,04:11 

U3U4 

60.2 

119.6 

2008-05-17,04:20 

U1U3 

101.5 

23.9 

2011-04-13,04:56 

U2U3 

44.1 

15.4 


2011-04-13,08:38 

U2U3 

46.2 

44.7 

2011-04-14, 09:35 

U2U3 

44.4 

48.9 

2011-04-16, 07:56 

U3U4 

62.3 

112.9 


2011-04-17,05:13 

U3U4 

53.7 

98.2 









HD 150193 

2007-06-03, 03:43 

U1U4 

128.0 

54.3 

2008-03-22, 05:24 

U2U4 

59.1 

43.1 

2008-03-23, 05:30 

U2U4 

60.9 

45.8 


2008-03-24, 07:36 

U1U4 

123.6 

48.3 

2008-03-24, 07:52 

U1U4 

125.4 

50.6 

2008-06-23, 06:21 

U2U3 

40.1 

51.5 


2013-04-29, 06:20 

U3U4 

62.1 

107.8 

2013-04-30, 05:06 

U1U4 

122.9 

47.5 

2013-04-30, 05:25 

U1U4 

125.1 

50.3 


2013-04-30, 06:04 

U1U4 

128.7 

55.5 

2013-05-02, 05:01 

U1U3 

100.6 

20.1 

2013-05-03, 04:20 

U1U2 

55.7 

9.0 


2014-04-16,06:50 

U3U4 

61.2 

105.7 









AK Sco 

2005-05-29, 03:27 

U1U2 

56.4 

13.6 

2005-05-29, 07:14 

U3U1 

91.8 

45.6 

2005-05-30, 02:48 

U3U4 

53.5 

91.1 


2005-05-30, 03:04 

U3U4 

55.2 

93.4 

2005-05-30, 08:44 

U3U4 

54.6 

151.6 





KK Oph 

2003-06-17,01:19 

U3U1 

100.3 

8.8 

2003-06-17,03:19 

U3U1 

102.2 

26.2 

2003-06-17, 03:40 

U3U1 

102.4 

28.8 


2003-06-17,05:55 

U3U1 

98.4 

41.4 









51 Oph 

2003-06-15,03:25 

U3U1 

101.3 

23.3 

2003-06-15,03:35 

U3U1 

101.5 

24.7 

2003-06-15, 06:54 

U3U1 

97.3 

42.4 


2003-06-15,08:23 

U3U1 

86.4 

44.5 

2003-06-16, 00:01 

U3U1 

98.9 

171.8 

2003-06-16, 00:09 

U3U1 

98.8 

173.2 


2003-06-16,02:17 

U3U1 

99.6 

14.1 

2006-05-16, 05:23 

U1U3 

101.3 

23.5 

2006-07-10, 03:55 

U2U4 

87.4 

86.7 

HD 163296 

2003-06-14, 03:13 

U3U1 

99.4 

17.7 

2009-07-07, 07:56 

HOGO 

16.8 

92.7 

2009-07-07, 07:59 

HOGO 

16.4 

93.1 


2009-08-14, 03:12 

EOGO 

14.2 

81.3 

2009-08-14, 03:53 

EOGO 

12.8 

84.4 

2009-08-15,00:16 

HOGO 

31.4 

66.1 


2009-08-15,00:31 

HOGO 

31.7 

67.8 

2009-08-15,02:26 

HOGO 

30.5 

78.1 

2009-08-15, 03:00 

HOGO 

28.9 

80.7 


2009-08-15,03:12 

HOGO 

28.2 

81.7 

2010-05-05, 05:03 

AOGl 

74.7 

104.2 

2010-05-05, 05:07 

AOGl 

75.5 

104.4 


2010-05-21,04:51 

HOEO 

43.7 

57.8 

2010-05-21,05:02 

HOEO 

44.5 

59.5 

2010-05-21,06:06 

HOH 

40.7 

146.3 


2010-05-21,07:05 

HOll 

40.6 

152.1 

2010-05-21,07:59 

HOll 

40.1 

158.7 

2010-05-21,08:52 

llEO 

60.6 

122.2 


2010-05-21,09:03 

llEO 

59.3 

124.1 









HD 169142 

2011-04-15,06:23 

U2U3 

45.1 

11.8 

2011-04-15,08:59 

U2U3 

46.6 

36.1 

2011-04-15,09:09 

U2U3 

46.6 

37.4 


2011-04-16, 08:49 

U3U4 

61.5 

106.0 

2011-04-17,07:52 

U3U4 

58.2 

100.2 

2011-04-17,08:00 

U3U4 

58.9 

101.1 


2012-06-06, 05:05 

U2U4 

87.4 

73.5 

2012-06-06, 05:13 

U2U4 

88.0 

74.6 

2012-06-06, 05:31 

U2U4 

88.8 

77.0 
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Table B.l. continued. 


name 

time 

(UTC) 

base 

Bp 

(m) 

P.A. 

(°) 

time 

(UTC) 

base 

Bp 

(m) 

P.A. 

(°) 

time 

(UTC) 

base 

Bp 

(m) 

P.A. 

(°) 

MWC 297 

2006-04-20, 07:55 

DOGO 

28.9 

68.8 

2006-04-20, 09:47 

DOGO 

32.0 

72.8 

2006-04-20, 10:29 

DOGO 

31.4 

73.2 


2006-04-20, 06:02 

DOGO 

20.1 

56.4 

2006-04-21,09:09 

DOGO 

31.7 

72.1 

2006-04-22, 09:59 

DOGO 

31.9 

73.1 


2006-05-23, 05:06 

AOGO 

52.6 

66.0 

2006-05-23, 06:52 

AOGO 

63.1 

71.8 

2006-05-25,05:12 

AOGO 

54.5 

67.1 


2006-05-25, 07:20 

AOGO 

64.0 

72.7 

2007-04-12, 07:52 

EOGO 

13.3 

66.3 

2007-04-12, 07:59 

EOGO 

13.5 

66.8 


2007-04-12, 08:22 

EOGO 

14.3 

68.5 

2007-04-17,06:50 

HOGO 

23.2 

61.9 

2007-04-18,08:36 

EOGO 

15.2 

70.5 


2007-04-19, 08:00 

EOHO 

43.3 

68.9 

2007-04-24, 06:54 

GIDO 

70.5 

129.5 

2007-04-24, 07:06 

GIDO 

70.9 

129.4 


2007-04-24,07:12 

GIDO 

71.1 

129.4 

2007-04-25,07:13 

HODO 

54.8 

67.3 

2007-04-25, 09:42 

GIDO 

66.1 

135.9 


2007-04-25, 09:54 

GIDO 

65.0 

137.1 

2007-05-08, 06:41 

GIDO 

71.5 

129.5 

2007-05-10, 05:45 

GlHO 

66.8 

172.9 

MWC 300 

2009-04-16,07:02 

EOHO 

36.5 

61.2 

2009-04-20, 05:56 

HOGO 

20.1 

52.7 

2009-05-03, 08:34 

GIDO 

69.2 

133.6 


2009-05-04, 06:21 

GlHO 

67.6 

174.3 

2009-05-04, 08:52 

GlHO 

68.4 

11.2 

2009-06-30, 05:24 

KOAO 

127.4 

73.2 


2009-06-30, 05:35 

KOAO 

126.8 

73.4 

2009-06-30, 06:06 

KOAO 

123.5 

73.8 

2009-07-01, 04:25 

AOGl 

90.0 

114.5 


2009-07-04, 05:56 

EOHO 

46.0 

73.8 

2009-07-04, 07:26 

EOHO 

37.7 

72.9 

2009-07-04, 07:34 

EOHO 

36.7 

72.6 


2009-07-05, 04:28 

EOGO 

16.0 

72.3 

2009-07-06, 05:57 

EOHO 

45.4 

73.9 

2009-07-06, 06:01 

EOHO 

45.2 

73.9 


2009-07-07, 05:01 

HOGO 

31.8 

73.3 

2009-07-07, 05:05 

HOGO 

31.7 

73.4 

2010-05-04, 09:34 

KOAO 

125.6 

73.6 


2010-05-04, 09:40 

KOAO 

124.9 

73.7 

2010-05-20, 05:49 

HOIl 

40.7 

142.8 

2012-05-01,07:34 

D0B2 

31.3 

21.0 


2012-05-02, 05:59 

B2A1 

10.2 

113.2 

2012-05-02, 09:55 

D0B2 

33.8 

33.3 





RCrA 

2004-07-09,09:10 

U3U2 

27.4 

63.5 

2004-07-28, 02:04 

U3U2 

46.6 

30.5 

2004-07-29, 23:49 

U3U2 

46.2 

9.4 


2004-07-30, 03:32 

U3U2 

45.6 

43.3 

2004-07-30, 05:33 

U3U2 

40.1 

55.5 

2004-07-30, 05:40 

U3U2 

39.6 

56.1 


2005-06-23, 07:33 

U1U4 

115.3 

76.1 

2005-06-25, 07:29 

U1U4 

114.6 

76.6 

2005-06-25, 07:39 

U1U4 

112.5 

77.8 


2005-06-26, 05:57 

U1U4 

127.5 

65.3 

2005-06-28, 04:40 

U3U4 

61.2 

105.0 

2005-06-28, 07:40 

U3U4 

58.8 

134.5 


2005-07-21,01:29 

U3U4 

53.3 

90.8 

2005-07-21,02:26 

U3U4 

58.7 

98.9 

2005-07-21,04:35 

U3U4 

62.3 

117.6 


2005-09-18,02:10 

U1U4 

111.2 

78.5 

2009-07-01,06:55 

AOGl 

88.3 

133.1 

2009-07-01,06:59 

AOGl 

88.1 

133.7 


2009-07-04, 08:36 

EOHO 

32.3 

106.1 

2009-07-04, 08:43 

EOHO 

31.4 

107.7 

2009-07-05,01:21 

EOGO 

14.0 

32.5 


2009-07-06, 04:43 

EOHO 

47.9 

70.9 

2009-07-06, 07:38 

EOHO 

37.8 

96.8 





TCrA 

2004-08-01,03:47 

U3U2 

45.0 

45.9 









HD 179218 

2003-06-16, 03:17 

U3U1 

59.7 

12.9 

2003-06-16, 03:31 

U3U1 

61.0 

16.6 

2003-06-16, 06:09 

U3U1 

84.7 

39.4 


2004-04-10, 07:43 

U3U2 

25.4 

27.0 

2004-04-10, 08:26 

U3U2 

29.1 

36.4 

2006-05-15,06:11 

U2U3 

29.1 

36.5 


2006-05-16, 08:01 

U1U3 

83.2 

38.8 

2006-05-16, 08:12 

U1U3 

85.0 

39.5 

2006-05-17, 05:39 

U3U4 

58.2 

118.0 


2006-06-11,06:42 

U3U4 

60.3 

105.7 

2006-06-14, 06:53 

U1U2 

48.2 

34.6 

2006-06-14, 08:20 

U1U2 

53.4 

36.8 


2006-07-09, 06:58 

U3U4 

44.4 

101.1 

2006-07-13,06:14 

U1U2 

52.8 

36.8 

2006-07-13,06:27 

U1U2 

53.5 

36.8 


2009-08-14,00:09 

EOGO 

11.6 

81.0 

2009-08-14, 00:47 

EOGO 

13.2 

79.7 

2009-08-14,01:20 

EOGO 

14.3 

78.3 


2009-08-14,01:51 

EOGO 

15.1 

76.9 

2009-08-14, 02:03 

EOGO 

15.3 

76.3 

2009-08-14, 02:36 

EOGO 

15.8 

74.5 


2009-08-14, 04:24 

EOGO 

15.4 

66.5 

2009-08-15,01:36 

HOGO 

29.7 

77.5 

2009-08-15,02:10 

HOGO 

31.1 

75.7 
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